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■ ABSTRACT 

o : 

i Context. A precise understanding of the relations between observable X-ray properties of galaxy clusters and cluster mass is a vital 

part of the application of X-ray galaxy cluster surveys to test cosmological models. An understanding of how these relations evolve 
with redshift is just emerging from a number of observational data sets. 

Aims. The current literature provides a diverse and inhomogeneous picture of scaling relation evolution. We attempt to transform 
these results and the data on recently discovered distant clusters into an updated and consistent framework, and provide an overall 
\q , view of scaling relation evolution from the combined data sets. 

I . Methods. We study in particular the most important scaling relations connecting X-ray luminosity, temperature, and cluster mass (M- 

T, L^-T, and M-Lx) combining 14 published data sets supplemented with recently published data of distant clusters and new results 
from follow-up observations of the XMM-Newton Distant Cluster Project (XDCP) that adds new leverage to efficiently constrain the 
scaling relations at high redshift. 

Results. We find that the evolution of the mass-temperature relation is consistent with the self-similar evolution prediction, while 
the evolution of X-ray luminosity for a given temperature and mass for a given X-ray luminosity is slower than predicted by simple 
(-H ' self-similar models. Our best fit results for the evolution factor E(zf are a = -1.04 ± 0.07 for the M-T relation, a = -0.23^^ for the 

O . L-T relation, and a = -0.93^ p for the M-L* relation. We also explore the influence of selection effects on scaling relations and find 

I that selection biases are the most likely reason for apparent inconsistencies between different published data sets. 

Conclusions. The new results provide the currently most robust calibration of high-redshift cluster mass estimates based on X-ray 
| luminosity and temperature and help us to improve the prediction of the number of clusters to be found in future galaxy cluster X-ray 

surveys, such as eROSITA. The comparison of evolution results with hydrodynamical cosmological simulations suggests that early 
preheating of the intracluster medium (ICM) provides the most suitable scenario to explain the observed evolution. 
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O ! 1 ■ Introduction (2008), Pratt et al. (2009), Mantz et al. (2009), Bohringeret al. 

. _ , , . , , (2010), Arnaudet al. (2009)), the understanding of the evolu- 

^ , Galaxy clusters have become important probes for the study tion of ^ ^ mds y faer redshifts . g kss ^ 

. ■ of the evolution of the large-scale structure and for the Resuhs ^ flow m that ide insi ht int0 the redshift 

ON _ test of cosmo ogical models (e.g. Borgam & Guzzo (2001 , fe d , = j feut ^ different duster ^ 

6t T f ■ ( °° I' ™^ 6t ?; • ( 2 S\ S ™ C 6t a ■ Provide inconsistent scenarios for the redshift evolution of the 

™ ' ^ lman 6t , S?' Jf- ( \ ™T 6t i scaling relations. Moreover, it is not trivial to compare the dif- 

(2009), Mantz et al. (2009), Bohringer et al (2010) and ferent 8 results since they have been partly deri ved for different 

provide an interesting laboratory to study galaxy evolution. cosmol ieSi for different denn itions of the scaling radius, and 

V i Nowadays the application of galaxy cluster surveys to cosmo- co - ed with different schemes of the self . similar scaling laws 

^ . logical studies is limited mainly by a lack of understanding of tQ if ^ evolution trend There f or e, this work makes 

& ■ galaxy duster properties and the precise scaling relations of an effort tQ ^ different ^ a uniform framework 

observables with cluster mass, in particular at higher redshifts. and combine ^ ayailable data tQ md ^ evoluti trend 

To date X-ray observations provide the most reliable and oyer ^ widest redshift baseline available 
detailed characterization of galaxy clusters and X-ray pa- 
rameters and are most widely used in cosmological galaxy The data sets used from the literature comprise the work 
cluster studies for several reasons: (i) X-ray luminosity is of Andersson et al. (2010), Pratt et al. (2009), Mantz et al. 
tightly correlated to the cluster mass (Reiprich & Bohringer ( 2009 )' Zhan 8 et aL ( 2008 ' 2007 )- Hlcks et aL ( 2008 )' 
(2002), Pratt et al. (2009)), (ii) bright X-ray emission is only Pacaudetal. (2007), Branchesi et al. (2007), O'Hara et al. 
observed for evolved clusters with a deep gravitational po- ( 2007 )' Vikhlimn et al. (2006), Maughanetal. (2006), 
tential well, and (iii) the X-ray emission is highly peaked, Arnaud et aL ( 2005 )' Kotov & Vikhlimn (2005), and Etton et al. 
minimizing projection effects. While great progress has been ( 2004 ) and we lnclude the data of an additional 15 clusters 
made in the characterizing clusters in X-rays at low redshifts from recent publications and the XMM-Newton Distant Cluster 
(Arnaud et al. (2005), Vikhlinin et al. (2006, 2009), Zhang et al. Pro J ect ( XDCP ' Fassbender (2008), Boehnnger et al. (2005)). 
The latter set of clusters significantly improve the statistics in 
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The derived evolution results are compared to the findings 
of hydrodynamical cosmological simulations assuming different 
heating and cooling scenarios and therefore allow an investiga- 
tion of the ICM thermal history. 

Tighter constraints on the evolution of the X-ray luminosity 
of clusters for a given mass also allow us to make refined predic- 
tions about the number of high redshift clusters to be observed 
in future X-ray surveys. We illustrate this in the context of the 
eROSITA mission (Predehl et al. 2010). 

The paper is structured as follows. In Sect. 2, we briefly de- 
scribe the cluster samples used and the way in which we have 
transformed these public results into a unified framework. We 
also outline the theoretical expectations of the scaling relations 
and an estimate of the selection bias inherent to our combined 
cluster sample. The local scaling relations and the results on their 
evolutionary trend with redshift are given in Sect. 3. In Sect. 5.1, 
these results are compared to the predictions of numerical simu- 
lation studies and the implications for the ICM heating scenario 
are discussed. In Sect. 5.2, we outline the impact of our results on 
the number of clusters to be detected with eROSITA and Sect. 6, 
we provide a summary and our conclusions. 

Throughout the article, we adopt a ACDM cosmology with 
(Q A , Q. M , H , w) = (0.7, 0.3, 70 km s" 1 Mpc" 1 , - 1). 



2. Data and data analysis 

2. 1 . The cluster sample 

The evolution of galaxy cluster X-ray scaling relations is in- 
vestigated by means of a combined cluster sample compiled 
from a number of recent publications. To enable constraints on 
the redshift evolution of scaling relations, the clusters were se- 
lected to cover a wide redshift range from local systems out to 
z - 1.46. We attempted to avoid the complications caused by the 
expected deviation of galaxy groups from the scaling laws for 
more massive clusters by applying an ICM temperature thresh- 
old of r m j n = 2 keV. 

Table 1 shows the publications used to compile the combined 
cluster sample. The acronyms listed there are used throughout 
the remainder of this paper when referring to the source publica- 
tions. The derived X-ray properties of clusters included in more 
than one subsample are compared in Appendix C. For these sys- 
tems, the results derived with an analysis scheme that is most 
similar to ours, hence requiring the least corrections are used for 
the combined sample. Selection biases and the intrinsic scatter 
in the cluster population about the mean relations ensure that a 
determination of scaling relation evolution in the redshift range 
up to z — 0.8 is challenging (see Sect. 2. 4). Therefore, an ex- 
tensive sample of high redshift clusters is crucial for this study. 
For this purpose, in addition to the subsamples listed in Table 1, 
high redshift clusters detected in the XDCP (Fassbender 2008), 
XMM-LSS (Pierre et al. 2004), XCS (Romeretal. 2002) sur- 
veys, the 2XMM catalogue (Watson et al. 2009), and by the 
South Pole Telescope (SPT) (AnderssonlO, Foley et al. (2011)) 
were included in the combined sample. 

Table A. 1 lists the galaxy clusters included in the combined 
sample consisting of 232 systems, of which 40 are at z > 0.8. 
The cluster temperature, X-ray luminosity, and mass used in the 
present analysis are listed in the table. The sample covers a range 
of cluster masses from about 5-10 13 M Q to 3-lO 15 M . For the dis- 
tant cluster sample, cluster masses derived by means of the Yx- 
M relation were not considered owing to their dependence on 
the Yx-M scaling relation and its redshift evolution. The z> 0.3- 
clusters used to place constraints on scaling-relation evolution 



cover an approximately uniform luminosity range over the entire 
redshift interval z = 0.3 - 1 .46 (see Appendix B) and by construc- 
tion, the distant cluster sample shows no strong morphological 
selection bias. 



2.2. Scaling theory 

Assuming that the evolution of the ICM during cluster forma- 
tion is governed solely by gravitational processes, clusters are 
expected to be self-similar objects whose X-ray properties and 
masses are connected by scaling relations predicted by the self- 
similar model (e.g. Kaiser (1986)). Since galaxy clusters have no 
clearly defined natural outer boundary, a fiducial radius within 
which cluster properties are considered has to be chosen. Scaling 
theory is used to define this radius in such a way that it de- 
scribes the same corresponding boundary for clusters of all sizes 
in the framework of the self-similar cluster structure model. 
In accordance with the homogeneous spherical collapse model 
of Gunn & Gott (1972) and detailed N-body simulations (e.g. 
Navarro et al. (1995)), the fiducial radius is defined to enclose a 
spherical region with a mean overdensity of A times the critical 
density of the Universe p C rit(z) 



3M(r < r A ) 
47rp crit (z)A ' 



(1) 



While this first-order self-similar model seems to describe the 
structure of dark matter haloes fairly well, additional gas physics 
including heating and cooling processes are needed to ex- 
plain the ICM structure and the resulting X-ray properties. 
Consequently, the local scaling relations have been found to 
differ from the self-similar predictions in some cases (see e.g. 
Pratt09), e.g. the Lx-T relation is steeper than expected. 

The self-similar model also predicts the evolution of scal- 
ing relations with redshift or lookback time. However, there are 
different schemes for defining the fiducial radius in order to en- 
close self-similar regions for clusters at different redshifts. One 
approach based on a spherical top-hat collapse model assumes 
that a galaxy cluster as it is observed has only recently formed 
at the given redshift and proposes adopting a redshift-dependent 
density contrast A- when defining fiducial radii, where A- can be 
expressed in terms of the density contrast at the virial radius at 
the cluster redshift 



A- = A(z = 0) 



A w ,-(z) 

A,„>(Z = 0)' 



(2) 



Bryan & Norman (1998) give an expression for A Vir (z) in a flat 
ACDM-cosmology of 



A vir (z) = 18;r 2 + 82[Q„,(z) - 1] - 39[Q„,(z) - if, 



(3) 



where ft m (z) = n M (l+z) 3 /£(z) 2 and E(z) = ^ . The expectation 
for the M-T, Lx-T, and M-Lx relation and their evolution with 
redshift in a model only taking into account gravitational effects 
is then 



Moc T 3/2 E( z y l A: 



1/2 
1/2 



L x oc T / -E(zW z 
MccLfzr 7 /^ 8 , 



(4) 
(5) 
(6) 



where Lx is the bolometric luminosity integrated out to the scale 
radius. This or related definitions of fiducial radii were used in 
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Publication 


Acronym 


Survey 


Instrument 


Cluster # 


z Range 


Andersson et al. (2010) 


Andersson 10 


SPT 


XMM/Chandra 


9 


0.4-1.1 


Pratt etal. (2009) 


Pratt09 


REXCESS 


XMM 


26 


< 0.2 


Mantz et al. (2009) 


Mantz09 


BCS, REFLEX 


Chandra/ROSAT 


42 


<0.3 


Zhang etal. (2008) 


Zhang08 


LoCuSS 


XMM 


37 


0.14-0.3 


Hicks et al. (2008) 


Hicks08 


RCS 


Chandra 


8 


0.6-0.9 


Pacaud et al. (2007) 


Pacaud07 


XMM-LSS 


XMM 


13 


< 1.05 


Branchesi et al. (2007) 


Branchesi07 


archival 


Chandra/XMM 


4 


0.25-0.46 


O'Haraetal. (2007) 


OHara07 


archival 


Chandra 


26 


0.29-0.82 


Zhang etal. (2007) 


Zhang07 


pilot LoCuSS 


XMM 


4 


0.27-0.3 


Vikhlinin et al. (2006) 


Vikhlinin06 


archival 


Chandra/ROSAT 


2 


< 0.2 


Maughan et al. (2006) 


Maughan06 


WARPS 


XMM/Chandra 


8 


0.6-1 


Arnaudetal. (2005) 


Arnaud05 


archival 


XMM 


5 


<0.15 


Kotov & Vikhlinin (2005) 


Kotov05 


archival 


XMM 


5 


0.4-0.7 


Ettori et al. (2004) 


Ettori04 


archival 


Chandra 


28 


0.4-1.3 



a number of recent publications on scaling relation evolution 
(e.g. Ettori et al. (2004) and Maughan et al. (2006)). The sec- 
ond commonly used definition is to measure cluster properties 
within regions with a redshift-independent value for the den- 
sity contrast A, which was applied e.g. by Pacaud et al. (2007) 
and Kotov & Vikhlinin (2005). The expected evolution of scal- 
ing relations in this framework is similar to Eqs. 4, 5, and 6, al- 
beit omitting the A bactors. We are currently unable to decide 
between the two approaches on the basis of observational data 
owing to the lack of sufficiently extensive and precise data sets 
needed to catch the subtle differences between the two models. 
In a forthcoming paper (Bohringer et al. 201 1), we explore this 
question by means of numerical simulations, finding that the 
recent formation approximation is imprecise and that the fixed 
overdensity approach (not including the A- factors) describes the 
simulation results more accurately than the formulae including 
this extra term. For this paper we decided to explore both ap- 
proaches and find that the differences are negligibly small for 
our conclusions. 



2.3. Homogenization scheme 

A number of corrections had to be applied to the subsamples in 
order to correct for the slightly different methods used by the 
authors. 

Throughout the cluster samples, we used both means of 
defining fiducial radii described in the previous section. Within 
these two schemes, various values of the mean overdensity are 
used. Most recent studies use a density contrast of 200, 500, 
or 2 500. The use of r2oo, which approximately corresponds to 
the virial radius, has the drawback that in many cases the area 
inside r2oo is not fully covered by the available X-ray data. In 
many clusters, r25oo typically corresponds to the most relaxed 
central part of the cluster and is also used in some publica- 
tions. The most common definition of the cluster radius that 
is also used in this work is (corresponding to either A- = 
500 • A vir (z) I A vir (z = 0) or a fixed A = 500). A large fraction 
of the clusters are well-relaxed within this radius, which is often 
well matched with the cluster region for which X-ray data are 
available. 

A rescaling scheme similar to the one used by 
Branchesi et al. (2007) was used to correct for the differ- 
ent definitions and values of density contrast. For this rescaling 
scheme, the cluster density profiles were assumed to follow the 
isothermal /?-model (Cavaliere & Fusco-Femiano 1976), since 



for the majority of the clusters included in the combined sample 
the shape parameters (3 and r c are known. In the framework of 
the ^-model, the mass enclosed within the radius r is given by 



M(r) = 



3(3k B T r(r/r c ) 2 
Gjjrn,, 1 + (r/r c ) 2 



(7) 



Using Eqs. 7 and 1, the radius corresponding to a density contrast 
of A, is given by 



Wk B T 



1 



fim p H(z) 2 A : 



For the /?-model, the fractional luminosity within r is 



L x (< r)/L tot = 



27^ o 
3/tot 



1 -(1 + (r/r c ) 2 ) 3/2 ~^ 



2/3-1 



(8) 



(9) 



where S o designates the central surface brightness and f tot the to- 
tal X-ray flux. Using these equations, a simple correction scheme 
can be applied to the cluster observables. First xi = r\ /r c , the ra- 
dius in units of r c for which the cluster properties are given is cal- 
culated using Equ. 8, where A- is calculated by means of Equ. 2. 
The radius corresponding to the density contrast to which the 
cluster properties will be rescaled, X2 = r2/r c , is calculated in the 
same way. Using Eqs. 7 and 9 we obtain the expressions 



M(x l )IM(x 2 ) 



L x (xi)/L x (x 2 ) 



\+x\) \l+x\ i 

1 - (1 + x]) 3 ' 2 ^ 
1 -(1 +*2)3/2-V 



(10) 

(11) 



The cluster properties are then multiplied by the correction fac- 
tors obtained by means of Eqs. 10 and 11. 

Some publications give cluster X-ray luminosities in either 
the 0.1-2.4 keV or 0.5-2 keV energy band. For the combined 
cluster sample, bolometric X-ray luminosities, that is cluster lu- 
minosities in the 0.01-100 keV-band, are used. The band lumi- 
nosities were converted to bolometric values by means of the 
X-ray spectral fitting package XSPEC (Mewe et al. 1985), as- 
suming a Mekal model with an ICM metallicity of 0.3Z Q . 

It has been shown (e.g. Markevitch (1998), Pratt09) that 
tighter scaling relations involving X-ray luminosity are obtained 
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by excising the cluster core region. For the combined sample 
selected in this work, however, the cluster observables for the 
local sample are given in a way that allows direct comparison 
to distant clusters for which core excision is not always feasi- 
ble. Therefore, luminosities in the entire r < f-500 aperture are 
considered, the emission from central regions is not excluded 
or replaced by any extrapolated profile to make the results of 
local samples comparable to the distant cluster studies. In the 
case of the ICM temperatures, obtaining a homogeneous dataset 
is less straightforward because in some studies of local clusters 
only core-excluded temperatures are given. Compiling a com- 
bined sample directly from studies of clusters where the core is 
included and others with core-excluded temperatures can lead 
to systematic bias. The source of this bias is the diversity of 
the central ICM temperature profiles that affects the difference 
between core-included and core-excised temperatures. In most 
cases, non-CC clusters have a flat central temperature profile but 
CC clusters with a rather pronounced cool core (CC) have a cen- 
tral temperature is between about one-third and one-half of the 
surrounding regions {e.g. Peterson et al. (2003)). The fraction of 
CC clusters therefore determines the magnitude of the induced 
bias. 

Subsamples with both core-included and core-excised tem- 
peratures available allow us to estimate the errors caused by the 
inhomogeneous measurement schemes. The sample of Pratt09 
consists of 3 1 low redshift clusters and was designed to be mor- 
phologically unbiased, hence its relative error introduced by us- 
ing core-excluded instead of core-included temperatures was 
found to be 7%. This value was added to the estimated temper- 
ature errors for samples with only core-excluded temperatures 
available (AnderssonlO, Zhang08, Zhang07, Arnaud05). For the 
cluster SPT-CL J2106-5844 (Foley et al. 2011), an emission- 
weighted core-included temperature of T = 8.5 keV was es- 
timated from the core-excluded temperature of 11.0 keV and 
the core temperature of 6.5 keV. As outlined above, the impor- 
tance of correct and homogeneous core exclusion throughout 
the combined sample depends on the abundance of CC clusters. 
For local systems, this abundance is found to be 40-70% (e.g. 
Hudson et al. (2010)). No consensus has emerged yet on how 
this CC fraction evolves with redshift. While Vikhlinin et al. 
(2009) find a decrease in the CC fraction from ~ 70% locally 
to ~ 15% at z > 0.5, Santos et al. (2008) find that the fraction 
of CC clusters at high redshift (z ~ 0.7 - 1 .2) is very similar to 
the local value (with an absence of very strong CCs). In sum- 
mary, using inhomogeneous temperature measurement schemes 
throughout the combined sample is believed to bias the evolution 
results only marginally. 

The values of most cluster properties depend on the assumed 
cosmological model. In all publications used to compile our 
combined cluster sample, the model of choice is the standard 
ACDM scenario. However, the value assumed for the Hubble 
constant Ho = 100/i km s _I Mpc~' varies slightly from h = 0.7 
to h = 0.73 within the list of source publications considered. 
The effects of this change slightly affect the values of the ba- 
sic cluster properties. To obtain comparable cluster subsamples, 
a common h of 0.7 is chosen and the necessary corrections are 
applied to the subsamples with a different h. 

Cluster mass depends on h as 

Mock' 1 , (12) 
while for the bolometric X-ray luminosities 

L x ccd 2 L och- 2 . (13) 



The cluster properties were rescaled accordingly, e.g. 
M70/M73 = 73/70 = 1.043 and L 7Q /L n = (73/70) 2 = 1.088. 

2.4. Selection bias estimate 

The cluster samples obtained in typical X-ray surveys are not 
strictly volume-limited but rather, at least approximately, flux- 
limited since the limited observation time, and detector area 
and sensitivity generally only enables the detection of objects 
brighter than a certain flux limit f m i n . Various selection biases 
complicate the analysis of these flux-limited samples and have 
to be taken into account when determining scaling relations and 
their evolution with redshift. For a well-controlled, homoge- 
neous survey, the cluster selection function, i.e. the probability 
of detecting a cluster with given properties taking into account 
the adopted observation strategy, can be modeled. A realistic se- 
lection function enables us to correct for selection effects in a 
more consistent and exact way than a survey for which only an 
approximate flux limit is known. An example of the application 
of this correction strategy can be found in Pacaud et al. (2007). 
However, the sample used in our work is highly heterogeneous, 
including clusters from numerous different surveys. For the ma- 
jority of these surveys, not all information necessary to recon- 
struct the survey selection function with high accuracy is avail- 
able. Therefore, a different strategy has to be applied to obtain at 
least a realistic estimate of the influence of selection effects for 
this sample. 

The approach adopted here consists of simulating a clus- 
ter population with as realistic as possible properties and se- 
lecting a cluster sample comparable to the observed one from 
this population. In this situation, in contrast to observed clus- 
ter samples, the characteristics of the underlying cluster popu- 
lation are known and the properties of the selected sample can 
be compared to those of the entire population to estimate se- 
lection effects. To probe the selection bias in the local scaling 
and the evolution of the L^-T relation, a temperature function 
nj{T, z) was assumed. Since no sufficiently exact measured tem- 
perature function was available, nj{T,z) was deduced from the 
more tightly constrained luminosity function «l(L, z) by multi- 
plying it with the determinant dL/dT and converting X-ray lumi- 
nosities to temperatures by means of the Lo.i-2.4keV-T relation. 
This method is only exact for the unrealistic case of no intrinsic 
scatter about the scaling relation, but provides a sufficiently ex- 
act approximation of the cluster temperature function to obtain 
a rough estimate of the selection bias. We note that for cluster 
discoveries and the related selection effects the luminosities in 
the X-ray observatory's detection band rather than the bolomet- 
ric X-ray luminosities used throughout this work are relevant. 
Since the majority of combined sample clusters originate from 
the ROSAT surveys, ROSAT band luminosities (0.1 - 2.4 keV 
observer frame) are used throughout this section. The number 
of clusters in a given redshift bin is therefore determined by the 
temperature function and the solid angle covered by the survey. 

Cluster luminosities were calculated assuming a non- 
evolving Lo.i-2.4keV-T relation since this evolution model rep- 
resents a fair first-order approximation to the observational 
data. The ROSAT band Lo.i-2.4keV-T relation given in Pratt09, 
£o.i-2.4keV = 0.078 • r[keV]) 2 - 24 1 44 ergs-', was used for this 
purpose. The luminosities were then displaced from the mean 
relation assuming a log-normal scatter of 0.25 dex (~ 60%), as 
suggested by the observed population and consistent with the 
value found by Pratt09. The chosen description is consistent with 
observations and seems reasonable when assuming a cluster's 
mass to be its basic property and taking into account the rather 
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Table 2. L^-T relations derived from the simulated local sam- 
ples. Sample 1 / m j n = 3 • 1CT 12 erg s" 1 crrT 2 , sample 2 / m i n 
I • 10" 12 erg s" 1 cm" 2 , sample 3 / m i n = 1 ■ 10" 13 erg s" 1 cm -2 . The 
input relation used to model the bolometric X-ray luminosities 
of the simulated cluster population is given for comparison. 



input relation 




L x = 


0.112 ■ (r[keV]) i ' ,, 10 44 erg s" 1 


sample 








-T relation 




1 


L x -- 


= (0.20 


± 0.03) • 


(r[keV])" 1±u 12 


10 44 erg s" 1 


2 




= (0.22 


± 0.05) • 


(r[keV]) 2 - 13±0 - 19 


lO^erg s" 1 


3 


L x -- 


= (0.09 


± 0.02) ■ 


(r[keV]) 2 - 78±02S 


10 44 erg s" 1 



tight correlation between mass and temperature compared to the 
greater intrinsic scatter in the L^-T relation. As for the observed 
cluster sample, a temperature threshold of T m i n = 2 keV was ap- 
plied to the simulated population. 

From this simulated cluster population, a flux-limited sample 
can be extracted including only clusters brighter than the limit- 
ing luminosity L m i n (z) resulting from the flux limit as 

L min (z) = 47Td 2 L f mi JK(z,T), (14) 

where di is the luminosity distance and K(z, T) the k-correction 
quantifying the relation between observer-frame band lumi- 
nosity and cluster rest-frame band luminosity, i.e. K(z, T) = 

^obs/ ^rest- 

The influence of selection effects on the observed local scal- 
ing relations was estimated by means of flux-limited samples se- 
lected from the simulated cluster population with limiting fluxes 
of /min = 3 • 10" 12 erg s" 1 cm" 2 (sample 1), 1 ■ 10" 12 erg s" 1 cm -2 
(sample 2), and 1 ■ 10" 13 erg s" 1 cm" 2 (sample 3), covering the 
redshift range used to fit the local relations, i.e. < z < 0.3. 
These flux limits approximately represent the range of limiting 
fluxes of the local cluster surveys used in this work. The survey 
areas were adjusted to provide sample sizes comparable to both 
each other and the local cluster sample used in this work, i.e. 
~ 100 clusters. The bolometric X-ray luminosities of the cluster 
population were assumed to follow the input L^-T relation of 
Equ. 17. We then fitted L^-T relations to the three samples us- 
ing the BCES(L|T) method (Akritas & Bershady 1996). Table 2 
summarizes the derived relations. 

In summary, the characteristics of the selection effects for lo- 
cal scaling-relation fits depend on whether the flux limit of the 
sample cuts away a significant fraction of the luminosity func- 
tion. For the faintest of the three samples described above (sam- 
ple 3), this is not the case, and consequently selection bias is 
negligible for this sample. Most of the clusters within the z<0.3 
sample used in this work were detected in surveys with rela- 
tively high flux limits, making the flux limit and not the applied 
temperature threshold the limiting factor at the low temperature 
and luminosity end of the cluster distribution. The bias in the 
observed sample is therefore expected to be comparable to the 
/min = 1 • 10 12 erg s" 1 cirr 2 and the / m i n = 3 • 10" 12 erg s" 1 cm" 2 - 
sample. As is clearly visible in the Lx-T relations fitted to these 
samples, the resulting bias for this range of flux limits is fairly 
insensitive to the exact limiting flux. This situation justifies a 
common bias estimate for the entire local sample used in this 
work. According to the simulated samples, the measured slope 
of the local L^-T relation is decreased only slightly by selec- 
tion bias, whereas the normalization is raised by almost 100%. 
The selection bias in the measured evolution of the Ly-T rela- 
tion with redshift displays similar trends, i.e. a higher normaliza- 
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Fig. 1. Evolution bias of the L^-T relation: Flux-limited sam- 
ples with /mm = 3 • 10" 12 erg s" 1 cm" 2 (green), / m i n = 1 ■ 
10" 13 ergs -1 cirr 2 (orange), /min = 1 TO" 14 erg s" 1 cm" 2 (black), 
and combined sample (blue curve). The bias for the combined 
sample rescaled to remove the bias in the local relation is plotted 
in red. The self-similar prediction for the evolution is plotted in 
grey. 

tion caused by selection effects, and is analyzed in greater detail 
throughout the remainder of this section. 

The simulated counterpart of the combined cluster sample 
used in this work was constructed from an all-sky survey with a 
flux limit of /min = 3 ■ 10" 12 erg s" 1 cm" 2 , representing the clus- 
ter samples based on the ROSAT All-Sky Survey (RASS), e.g. 
the REFLEX (Bohringer et al. 2004) survey. A second sample 
with a flux limit of f m \ n = 1 ■ 10" 13 erg s" 1 cm" 2 and an area of 
400 square degrees was added, representing the clusters from 
ROSAT PSPC -based surveys such as WARPS (Perlman et al. 
2002) or the 400sd (Burenin et al. 2007) survey. As a third com- 
ponent, a sample with / m j n = 1 ■ 10" 14 erg s" 1 cm" 2 , an area of 80 
square degrees and a minimum redshift of z m in = 0.8 was added, 
corresponding to current serendipitous surveys such as XDCP 
(Fassbender 2008). 

With this simulated sample at hand, the influence of bias on 
the evolution of the L^-T relation can be estimated. To achieve 
this, we calculated the logarithmic mean of the bolometric 
cluster luminosity divided by the luminosity resulting from the 
bolometric L^-T relation at the cluster temperature in redshift 
bins. This number quantifies the selection bias, with a value of 
1 corresponding to no selection effects. The value of the bias 
curve in a redshift bin with a width of Az and centered on the 
redshift z is therefore given by 

^clusters lo e( jfz) 

B( z )= 10 , (15) 

where the subscript "clusters" designates all clusters within [z - 
|Az; z + |Az] included in the flux-limited sample and A and b 
quantify the normalization and slope of the local Lx-T relation 
and were set according to Equ. 19. 

The bias curve deduced for the simulated combined sample 
is indicated with a blue-dashed line in Fig. 1. However, in the 
redshift range used to fit "local" scaling relations, 0.05 < z < 0.3, 
selection effects are already non-negligible, visible in a clear de- 
viation of the bias curve from 1 . A value of 1 on the vertical axis, 
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theoretically corresponding to no bias, therefore already includes 
the bias present in the local sample used to fit the Lx-T relation. 
This effect is redshift-independent because the same Lx-T rela- 
tion is used for clusters at all redshifts to compare their measured 
luminosity to the expected one. To distinguish the component of 
the selection bias that may mimic evolution of the Lx-T rela- 
tion, the effects of the local bias have to be taken into account. 
The local scaling relation bias leads to a division of cluster lumi- 
nosities by a higher expectation value than the unbiased relation. 
To determine the evolution bias relative to this higher expecta- 
tion value and not relative to the underlying cluster population, 
the bias curve has to be divided by the mean bias in the redshift 
range that was used to fit the "local" scaling relation. The bias 
curve was therefore rescaled by a factor of 0.8, corresponding to 
the red-dashed curve in Fig. 1. To summarize, the rescaled bias 
curve shows the additional redshift-dependent bias of the loga- 
rithmic mean of the bolometric cluster luminosity divided by the 
luminosity resulting from the local Lx-T relation relative to the 
already bias-affected local sample. 

The bias curves of Fig. 1 have various characteristics that can 
easily be explained in terms of the underlying cluster sample. 
The non-rescaled bias is generally greater than 1 because for 
any flux-limited sample in the presence of scatter more clusters 
below the mean relation than above it are too faint to be included 
in the sample. With increasing redshift, the fraction of the cluster 
population that is excluded for being too faint increases, causing 
an increase in the fraction of clusters above the mean Lx-T rela- 
tion that have no counterpart below the relation. Hence, for a sin- 
gle flux limit the bias increases with redshift. This trend is visible 
in the bias curves for the three flux-limited subsamples, / m ; n = 3 ■ 
10~ 12 erg s cm _2 (green), / m i n = 1 • 10~ 13 erg s _1 cirT 2 (orange), 
and / m i„ = 1 • 10~ 14 erg s _1 cm~ 2 (black). Up to a certain red- 
shift threshold, the influence of selection bias on the subsam- 
ples is negligible because no significant part of the cluster pop- 
ulation is excluded from the sample owing to the flux limit. 
Naturally, the unbiased redshift range increases with the sur- 
vey sensitivity from about z ~ 0.05 for the sample with /min = 
3 ■ 10~ 12 erg s _1 cirT 2 to z ~ 0.2 for the sample with / m j n = 
1 ■ 10~ 13 erg s _1 cirT 2 and z~ 0.6 if f m i„ = 1 • 10~ 14 erg s _1 cm -2 . 
Surveys with a high limiting flux display a faster increase in bias 
with redshift than deeper surveys. 

The characteristics of the mean bias curve for the combined 
cluster sample (blue curve in Fig. 1) are affected by the domi- 
nant contribution in terms of cluster detections with increasing 
redshift shifting from the all-sky / m j n = 3 • 10~ 12 erg s _1 cm -2 - 
survey to the more sensitive but smaller solid-angle serendipi- 
tous surveys. The combined bias curve therefore decreases when 
the contribution from the f m [ n = 1 • 10~ 13 erg s _1 cirT 2 -survey be- 
comes dominant at about z ~ 0.15. At z ~ 0.8, the curve drops 
again because from there toward higher redshift the contribu- 
tion of the / mm = 1 • 10~ 14 erg s _1 cirT 2 -surveys dominates. From 
z = 0.8 toward higher redshift, the bias increases again as out- 
lined above. 

The bias in the M-Lx relation was not determined indepen- 
dently but based on the results of the Lx-T relation. Cluster 
masses were set according to the local M-T relation (Equ. 16) 
and the ICM temperatures of the simulated cluster sample out- 
lined before. The logarithmic mean of the cluster masses divided 
by the masses expected from the M-Lx relation (Equ. 18) was 
then calculated for the simulated sample to obtain a bias curve 
analogous to the one derived for the Lx-T relation. By construc- 
tion, the bias curves for the M-Lx relation show the same fea- 
tures as those for the Lx-T relation. However, the bias curve is 
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Fig. 2. Redshift-dependent bias in the measured intrinsic scatter. 
Red-dashed line: true intrinsic scatter of the simulated samples 
Black: fitted intrinsic scatter for / mm = 1 
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inverted and the selection effects generally lead to an underesti- 
mation of the mean mass for a given luminosity. 

The intrinsic scatter in the cluster observables about the 
mean scaling relations is of great importance to the interpreta- 
tion of results about scaling relation evolution. While the intrin- 
sic scatter for local clusters is roughly known (see Pratt09), its 
redshift-dependence has not been well-constrained. Measuring 
the intrinsic scatter at high redshifts is challenging since in addi- 
tion to the intrinsic scatter, the observed scatter in the evolution 
plots of Sect. 3.2 has a number of other causes. First of all, mea- 
surement errors naturally have an influence on the observed scat- 
ter. Furthermore, increased scatter can also result from a change 
in the scaling-relation slope with redshift. Finally, selection ef- 
fects of flux-limited cluster samples may have an influence on 
the observed scatter. 

To estimate this bias in the Lx-T relation, we measured the 
intrinsic scatter of various simulated cluster samples. The sim- 
ulated cluster population has an intrinsic scatter in luminosity 
of 0.25 dex (~60%). From this population, flux-limited samples 



with/ min = M0 erg s" 



and/ min = M0 14 erg s" 



were selected. The intrinsic scatter of the samples was then fit- 
ted in redshift bins centered around z- 0.1, 0.3, 0.7, 0.9, 1.1, and 
1.3. The results are shown in Fig 2. 

For all subsamples, the fitting routine generally underesti- 
mates the real intrinsic scatter. The fitted scatter for local clus- 
ter samples is about 0.22 dex and comparable for the two sam- 
ples. As more and more luminous clusters are excluded by the 
flux limit at higher redshifts, the fitted scatter shows a decreas- 
ing trend. As expected, the decrease is more rapid with red- 
shift for the / m j n = 1 • 10~ 13 erg s _1 crrT 2 -sample than the deeper 
/min = 1 • 10~ 14 erg s crrT 2 -sample. For the most distant sub- 
samples at both flux limits (z~ 1.3), the estimated scatter is 0.12 
and 0.14 dex, respectively, i.e. the selection bias causes the scat- 
ter to be underestimated by about 50%. 
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3. Results 

3.1. Local scaling relations 

The results for scaling relations fitted to observed clus- 
ter samples may differ significantly depending on the fitting 
scheme used. In this work, we use the BCES fitting method 
(Akritas & Bershady 1996) that has been widely used in re- 
cent studies and correctly accounts for intrinsic scatter about 
the mean relation and inhomogeneous measurement errors. 
However, several slightly different variations of this method ex- 
ist. When choosing one of these alternatives, it is most important 
to distinguish between a fundamental independent and depen- 
dent variable. For the cosmological applications related to the re- 
sults of this study, cluster mass is the fundamental property, and 
luminosity and temperature take the role of dependent variables. 
Therefore for the M-T and M-L x relation, the BCES(T|M) and 
BCES(L|M) method is used in this work. Owing to the large 
intrinsic scatter in the M-L x relation, the ICM temperature dis- 
plays a tighter correlation with cluster mass than the X-ray lumi- 
nosity. As a consequence, for the L x -T relation the BCES(L|T) 
scheme is used. 

Only clusters with z < 0.3 were considered for the local fit. 
This redshift threshold was chosen to be large enough to include 
a sufficiently large number of systems and improve the quality of 
our statistical analysis but small enough to keep evolution effects 
to a negligible level. To determine the influence of evolutionary 
effects, the fits were repeated with a maximum redshift of z = 0.2 
instead of z — 0.3. The results derived with this lower redshift 
threshold are fully consistent with the sample at z < 0.3, that is 
even though the expected evolution factor is non-negligible at 
z = 0.3, no obvious evolutionary effects on the scaling fit for the 
entire local sample are observed. However, the statistical errors 
increase significantly with redshift because of the smaller sample 
size. The cluster properties were not rescaled by any assumed 
evolutionary model before the fit. 

The scaling relations derived from the combined local cluster 
sample are 

M = (0.236 + 0.031) • (r[keV]) L76±OO8 lO 14 M , (16) 
L x = (0.112 + 0.031) • (r[keV]) 2 53±015 1 44 erg s" 1 , (17) 



M = (1.191 + 0.104) ■ (L x [10 44 erg s -1 ])' 



6±0.04jq14^ 



(18) 



where L x is the bolometric X-ray luminosity and all properties 
are considered to be within rsoo. However, for the analysis of 
scaling relation evolution with redshift, the relations derived by 
Pratt09 were considered instead of the fitted scaling relations be- 
cause the former were derived from a more homogeneous sam- 
ple with well-known selection criteria, that have smaller relative 
errors. The L x -T relation of Pratt09 



L x = (0.079 + 0.008) • (r[keV]) 2 - 70±a24 1 44 erg s" 1 



(19) 



is consistent with our result in terms of both slope and normal- 
ization. No M-T relation is provided in Pratt09, but a fit to their 
sample with the BCES(T|M) method yields 

M = (0.291 + 0.031) ■ (r[keV]) L62±0 08 1 14 M o , (20) 

which is consistent with our result. Pratt09 provide an L-M in- 
stead of an M-L x relation. Fitting the inverse relation to their 
sample leads to 

M = (1.39 + 0.07) • (L x [10 44 erg s" 1 ]) - 5 ^ 003 ^ 1 ^, (21) 



which differs slightly from our result at the <2<x level. 

The local scaling relations derived by means of the differ- 
ent fitting methods and the relations by Pratt09 are shown in 
Appendix D. 

3.2. Evolution constraints from the combined cluster sample 

The central goal of this work is to obtain a clearer understanding 
of the redshift evolution of X-ray scaling relations. The figures in 
this section help us to visualize these evolutionary trends. They 
show the redshift-dependent distribution of cluster properties di- 
vided by the expected value assuming local scaling relations. For 
the Lx-T relation, this means that all cluster luminosities for in- 
stance are divided by L z= o(T), that is the luminosities inferred 
from the local L x -T relation at the measured cluster tempera- 
ture, and that we plot the quantity L ^ r) as a function of redshift. 

Plotting the data this way, the properties of the local clusters 
have an approximately log-normal scatter around 1 . A change in 
the normalization of the scaling relation with redshift translates 
into similar scatter around a different mean value in the evolu- 
tion plot, whereas a change in the slope would result in a larger 
scatter around the mean value. We note, however, that this is not 
a very suitable test for changes in slope as a greater intrinsic 
scatter at earlier times also leads to larger scatter in the evolu- 
tion plot and as outlined in Sect.2.4, selection biases may lead 
to an underestimation of the scatter of up to 50%. To investigate 
changes in the slope of high-z scaling laws, relations were fitted 
to the available z > 0.8-clusters by means of the BCES method. 
Owing to the small sample size, however, the errors in the result- 
ing relations are too large to allow independent constraints. 

The self-similar model of cluster formation predicts the 
slopes of X-ray scaling relations to be redshift-independent and 
the normalization to vary in proportion to powers of E(z) in the 
case of fixed overdensity. To test these predictions, a power-law 
E(z) a was fitted to all cluster data points with z> 0.3 for the M-T 
relation, i.e. the unknown exponent a in 



M ohs 
M z= o(T) 



= E(zf 



(22) 



was constrained by fitting M M ^ s r) versus E{z) in log-log-space. 
In the redshift range 0.3 < z < 0.6, an estimate of the influence 
of selection biases on the evolution results is challenging, since, 
in contrast to the more distant systems, the properties of these 
systems were obtained from both shallow surveys with a very 
high influence of selection biases but also deeper recent surveys. 
Therefore, this redshift range was excluded from the evolution 
fits for the L x -T and M-L x relation, for which selection biases 
play a more critical role than for the M-T relation. The excluded 
redshift range contains 45 sample clusters, while 65 clusters at 
z > 0.6 are used for the evolution fit. Owing to the large number 
of luminous clusters from shallow surveys, whose lack of depth 
introduces the largest bias into the overall sample, an inclusion 
of this redshift range in the evolution fit would lead to a positive 
evolution result that does not trace the observed evolution for 
more distant clusters from deeper surveys. To account for the 
variations in the number density of clusters with redshift in our 
sample, i.e. to avoid the result being exclusively determined by 
the large number of relatively low-redshift clusters with small 
errors, the data points were weighted by the inverse number 
of clusters in the corresponding redshift bin (Az = 0.1). Our 
evolution results are summarized in Table 3. 
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Fig. 3. Redshift evolution of the M-T relation. Black-dashed line: self-similar prediction (oc E(z) '). Continuous red line: best-fit 
evolution (oc E(zT 1M±0 07 ). 



Table 3. Evolution results based on the combined cluster sample. 
First column: Scaling relation. Second column: Observed evolu- 
tion of scaling relations, bias effects have been accounted for by 
greater uncertainties. Third column: Evolution results including 
a tentative selection-bias correction. Fourth column: Self-similar 
expectations and self-similar predictions 



relation 


observed evolution 


bias-corrected 


self-similar 


M-T 


oc E( Z r iM±m 




oc E(zY l 


Lx-T 


OC £( Z )-°' 23 -0.62 


£ (z) -0.6 5± 0.13 


oc E(z) +[ 


M-L x 


OC E(l) -0-12 


E(z) -om±QA2 


oc E(z)- 114 



We first discuss the M-T relation, which is expected among 
all relations to most closely follow the self-similar predictions. 
Figure 3 shows the redshift evolution of the M-T relation for the 
combined cluster sample. The best fit to the data corresponds 
to a redshift dependence of the normalization proportional to 
E(z)~ lM±om , which is consistent with the self-similar prediction 
of E(z) . We note that selection bias is not taken into account 
in this plot. However, for the the M-T relation this is not as im- 
portant as for scaling relations including cluster luminosity. An 
analysis similar to the one shown in Fig. 3 was performed with 
cluster radii defined with a variable density contrast A, instead 



of a fixed A at cluster redshift. The results are consistent within 
the errors with those of Fig. 3, i.e. using the redshift-dependent 
density contrast does not significantly influence our results about 
the evolution of the M-T relation. The M-T relation fitted to the 
z > 0.8-clusters with sufficiently good X-ray data has a slope of 
Mocr L59±0 - 45 , which is fully consistent with the local result. 
Figure 4 shows the redshift evolution of the L^-T relation for the 
combined cluster sample. The best-fit relation for the evolution 
is E(zy°' 23±0A2 , that is there is a slightly negative evolution. This 
result is clearly inconsistent with the self-similar prediction that 
the normalization increases with redshift in proportion to E(z) +1 . 

The evolution result was uncorrected for the estimated selec- 
tion bias (see Sect. 2. 4) because this bias estimate relies on a toy 
model that is only approximately comparable to the real cluster 
sample. The error budget was instead increased by the estimated 
bias, i.e. the confidence region was enlarged by the size of the 
estimated bias (see Fig. 1) in the direction of the supposed bias 
correction. This led to a final evolution result of E(z) -°- 62 . As 
expected, applying an approximate bias correction based on the 
rescaled bias curve of Fig. 4 before the fit as a test of the influ- 
ence of selection biases results in an even more negative evolu- 
tion result of £( z )-0-65±o.i3 

The slope of the L^-T relation fitted to the z > 0.8-clusters 
with sufficiently accurate X-ray data available is L oc 7- 3 - 12± 0-37 
This result is slightly steeper but still consistent within the errors 
with the local slope derived by Pratt09. Owing to the small clus- 
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Fig. 4. Redshift evolution of the Lx-T relation. Black-dashed line: self-similar prediction (oc E(zJ). Continuous red line: best-fit 
evolution (oc E(z)~ 02i -<>^). Red-dashed line: estimated mean bias for the combined sample rescaled to remove the effects of bias in 
the local scaling relation. 



ter sample and the large errors, this result heavily depends on 
the fitting method used and therefore provides no significant ev- 
idence of a steepening of the high redshift Lx-T relation. As for 
the M-T relation, the use of a redshift-dependent density con- 
trast A z instead of a fixed A leads to comparable results with 
similar scatter about the mean relation. 

Figure 5 shows the redshift evolution of the M-Lx rela- 
tion for the combined cluster sample. The best-fit relation is 
E(z) _0 - 93±012 , i.e. negative evolution as predicted by the self- 
similar model. However, our result is significantly less steep 
than the self-similar prediction of E(zy 7 ^ 4 . As for the M- 
T relation, the estimated selection bias is taken into account 
in the error budget and leads to the evolution being propor- 
tional to oc E(z) -0 - 93 -o'i2. This observed evolutionary trend is 
close to the one expected after combining the results for the 
evolution of the Lx-T and M-T relations, which would be 
oc E(zy 90 . Applying an approximate bias correction to test the 
influence of selection biases as for the Lx-T relation results in 
a slightly less negative evolution fit of E(z)~° M:k012 . \ n Pratt09, 
a bias-corrected local Lx-M relation is provided. Using the in- 
verted bias-corrected BCES(L|M)-relation (M = (1.64 + 0.07) ■ 
(Lx[10 44 erg s _1 ]) 0t52±003 ) and correcting for the total estimated 
evolution bias (not the curve rescaled to remove the effects of 
bias in the local scaling relation) leads to an evolution result of 
£(z) _0 - 90 -oi5 . The estimated errors in this result include statisti- 



cal errors and an estimate of the systematical error caused by an 
inexact bias correction. 

The slope of the fitted high redshift M-Lx relation is M oc 
L°- 70±0 - 21 , which is consistent with the local result. Using the A z - 
scheme instead of a redshift-independent density contrast again 
leads to similar results. The observed scatter in cluster properties 
about the mean relation for the high redshift clusters is consis- 
tent with the local scatter in all three relations. The real scatter 
about the Lx-T and M-Ly relation for distant clusters may be 
up to a factor of two larger than the observed result because of 
the influence of selection biases (see Sect. 2.4). However, owing 
to conservative error estimates, no constraints on the intrinsic 
scatter in cluster properties can be placed based on the measured 
total scatter for the distant cluster sample. 

4. Discussion 

4. 1 . Stability of results 

The results on scaling relation evolution presented in the previ- 
ous section were obtained by means of a number of input as- 
sumptions and results of preceding studies that have an influ- 
ence on the obtained results and may introduce additional errors. 
Throughout the remainder of this section, we briefly discuss the 
stability of the results under these assumptions. 
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Fig. 5. Redshift evolution of the M-Lx relation. Black-dashed line: self-similar prediction (oc E{z) 7 ^ 4 ). Continuous red line: best-fit 
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The assumed local scaling relations have a direct influence 
on the observed evolution. For our analysis in Sect. 3.2, the red- 
shift evolution of scaling relations was determined using the lo- 
cal scaling relations of Pratt09. Using the relations derived for 
the entire local combined sample (see Sect. 3.1) instead of these 
results does not lead to fundamentally different findings on the 
evolution coefficient. For the M-T relation, using the relation 
derived from the combined sample leads to a best-fit evolution 
of E(z)~ l 15±0 06 , which is only slightly different from the result 
for the Pratt09 relation (E(zy lM±om ). Using the L x -T relation 
fitted to our sample implies an evolution result of £ , ( z )-°- 36±0 - 12 
instead of E{zT 02 ^ and for the M-Lx relation £( z )-i06±o.i3 
instead of £ , (z) _0 - 93 -o i2, both of which are fully consistent with 
the results presented above for both relations. 

An incomplete or incorrect homogenization scheme applied 
to the different subsamples naturally influences the evolution re- 
sults. However, the combined cluster sample provides no hints 
that this might be a major problem (see Appendix C). Incorrect 
homogenization would be visible as larger scatter about the 
mean behavior in the cluster sample or different evolutionary 
trends for different subsamples. Taking into account selection 
biases, these significant trends are not observed for the cluster 
sample (see Sect.4.2). 

Our study again highlights the importance of selection bi- 
ases when investigating scaling relation evolution and the prob- 



lems inherent to small cluster samples over a limited redshift 
range. Although the simulated cluster sample used to estimate 
bias effects in this study is only a rough approximation of the 
true situation, it reveals the apparent evolutionary trends caused 
by selection effects, which have been taken into account in the 
estimated errors. Despite the lack of knowledge about the exact 
effects of selection bias in a highly inhomogeneous combined 
cluster sample such as ours, at least a fair estimate of the influ- 
ence on the evolution results can be given. For both the M-Lx 
and the Lx-T relations, a bias correction of the evolution results 
would render the difference to the self-similar predictions even 
more significant. Our finding about the inconsistency with the 
self-similar model can therefore not be attributed to selection ef- 
fects. 

The X-ray properties of the systems within a cluster sam- 
ple and their evolution with redshift might depend on the clus- 
ter selection strategy. While the subsamples of clusters selected 
by means of the Sunyaev-Zel'dovich effect (SZ) and their opti- 
cal/infrared (IR) properties are not sufficiently extensive to place 
independent constraints on scaling relation evolution, the X-ray, 
SZ,and optical/IR-selected subsamples display no obvious dif- 
ferences in their evolutionary trends. 

In recent studies, we note that different definitions of the den- 
sity contrast A and the resulting cluster radii have been used. 
However, the choice of either a fixed A at the cluster redshift or 
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a redshift-dependent A, has no significant influence on the deter- 
mined evolution results. 

4.2. Comparison to published results 

We now briefly discuss previously published results on scaling 
relation evolution comparing these to our findings: 

Kotov05: For the M-T relation, their evolution result of oc 
E(z)~ 0&i±023 is consistent with both our results and the self- 
similar expectation. For the L^-T relation, Kotov05 find the 
normalization to be oc (1 + z) L8±a3 , i.e. a positive evolu- 
tion that is even steeper than self-similar. This trend is eas- 
ily visible for their sample in Fig. 4. The seven clusters of 
the Kotov05 sample cover a redshift range of 0.4 < z < 0.7. 
Comparing this subsample with the remainder of the com- 
bined sample in this redshift range reveals that instead of 
describing the true evolution, this result can be attributed to 
selection effects, in addition to both limited sample size and 
redshift range. We note that the evolution result derived in 
Kotov05 approximately traces the bias curve for the sample 
with / m j n = 1 ■ 10~ 13 erg s _1 crrT 2 visible in Fig. 1. Although 
no uniform flux limit can be assigned to this archival sam- 
ple, /min = 1 • 10" 13 erg s _1 cirT 2 represents an appropriate 
estimate of the mean flux limit, i.e. the bias estimate sug- 
gests that after correcting for selection effects the observed 
evolution should be close to zero. 

Maughan06: Their result for the L^-T evolution is oc (1 + 
z)° 8±0 4 when using the local relation of Arnaud & Evrard 
(1999) and oc (1 + z ) 7±0 4 for the Markevitch (1998) re- 
lation, i.e. a slightly positive evolution. Using the scaling 
relation derived by Pratt09 instead removes this evolution- 
ary trend, causing the evolution of the Maughan06 subsam- 
ple (oc E(zy 0M±Q31 ) to be consistent with the results of our 
study. 

OHara07: Their result about the evolution of the core-included 
Lx-T relation within rsoo is E(z) 1 (l +z) _1 ' 90 -i i.e. consistent 
within the errors with our result. 

Branchesi07: No significant evolution of the L^-T relation in 
the redshift range 0.3 <z< 1.3 was observed. However, the 
normalization for this redshift range was found to be about 
a factor of ~ 2 higher than suggested by the local relations 
of Markevitch (1998) and Arnaud & Evrard (1999). The dif- 
ferences from the results derived in this work can mostly be 
attributed to the use of a different local scaling relation. 

Pacaud07: Since a detailed selection function was derived in 
this work, this analysis provides an important insight into 
the influence of selection effects. Before correcting for those, 
their result about the evolution of the Lx-T relation is 
roughly similar to that of Kotov05. Afterwards, they obtain 
an evolution factor of oc E{z) l {\ + z)~ om -°'- 55 , i.e. slightly less 
than self-similar. This result is still marginally inconsistent 
with ours. We note, however, that a significant part of their 
sample was not included here because the temperatures were 
below 2 keV. 

Ettori04: Their inferred evolution of the L^-T relation varies 
drastically depending on whether clusters below z < 0.6 are 
included in the fit. Using the local relation of Markevitch 
(1998), they find the evolution to be oc (1 + z f-(>2±o.2S for the 
entire sample and oc (1 + z ) O O4 ±0-33 jf on jy c i us t ers with z> 0.6 
are considered. Using the relation of Pratt09 leads to results 
consistent with those from the entire combined sample at 
least for the z > 0.6-clusters (oc £(z)-°- 60±019 ). The marked 
difference between the clusters at z < 0.6 and the more dis- 



tant systems can probably be attributed to selection effects, 
that is to the fact that the low-z sample consisting mostly of 
archival clusters detected in relatively shallow observations, 
as discussed above for the results of Kotov05, in accordance 
with which Ettori04 find the M-T evolution to be consistent 
with the self-similar prediction. 
Vikhlinin et al. (2009): When applying X-ray galaxy cluster 
data to cosmological tests these authors also provide a new 
evaluation of the L^-M relation and its evolution with red- 
shift in their Equ. 22. In contrast to our calculations, this re- 
lation is determined for luminosities in the 0.5 - 2 keV band. 
Using the results of Pratt about the difference between the 
bolometric and 0.5 - 2 keV band scaling relations for com- 
parison, we find that the results of Vikhlinin et al. and ours 
are in very good agreement for the zero redshift relation. 
Analyzing the redshift-dependent term in the relation infer- 
ring Lx from M, we find a term of £ , ( z )i- 85±0 42 f or Vikhlinin 

i 70+041 

et al. and a term of E(z) by inverting our Equ. 26. 

There is again good agreement within the large error bars. 
Taking the mean trends of both relations, we find a differ- 
ence in the evolution parameter of ~ 5% at redshift z — 0.8 
and ~ 10% at z = 1.5. Therefore, the difference between the 
two relations is smaller than can be detected with any current 
data set. 

Leauthaud et al. (2010): In their study, the M200-LX relation 
for 206 X-ray-selected galaxy groups from the COSMOS 
survey is determined in the redshift range 0.22 < z < 0.9 by 
means of stacked weak-lensing mass estimates. The derived 
evolution of oc E(z)~ inl 9 is consistent with our result. 

In summary, the present study provides a clearer picture of 
the scaling-relation evolution in three aspects. The first is that 
the use of the very early results by Markevitch (1998) and 
Arnaud & Evrard (1999) as a local reference of the scaling rela- 
tions introduces a positive evolutionary trend as already shown 
in Branchesi07. Using more recent results with higher quality 
statistics, and in particular the use of the weakly biased (not flux- 
limited) REXCESS sample results by Pratt09 removes most of 
this trend. The second aspect is that the overview of a larger set 
of cluster samples with different biases leads to the identification 
of bias effects from high flux-limits. Thirdly, the extension of the 
redshift range to newly detected high-redshift clusters increases 
the leverage for the evaluation of evolutionary effects. 

5. Implications of results 

5. 1. Implications for the thermal history of the ICM 

The observed modification of the evolution of scaling relations 
with respect to the self-similar model is caused by changes in 
the thermodynamic al state of the ICM. Therefore, different evo- 
lutionary trends are the signature of different histories of heating 
and cooling processes in the clusters. Comparing observational 
results with the predictions made by hydrodynamical cosmolog- 
ical simulations permits us to assess whether the heating scheme 
implemented in the simulations is realistic. A variety of exist- 
ing simulations taking into account different sources of heating 
and assumptions on the time evolution of non-gravitational heat- 
ing and cooling allow a constraint of the most realistic heating 
scheme. 

Previous attempts these comparison analyses were compli- 
cated by the inconsistent observational results about scaling re- 
lation evolution. Owing to updated local scaling relations based 
on morphologically unbiased cluster samples and the availabil- 
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ity of a larger cluster sample over a wider redshift range, evolu- 
tion constraints derived from the combined cluster sample used 
in this work permit a meaningful comparison. The simulations 
considered on that account are the Millennium Gas Simulations 
(MGS, Staneketal. (2010), Short et al. (2010)), a series of 
hydrodynamical resimulations of the original dark-matter-only 
Millennium simulations (Springel et al. 2005). The MGS runs 
incorporate a (500/i _l Mpc) 3 -volume with the same initial con- 
ditions and cosmological model as the original Millennium run, 
but include an equal number of gas particles in addition to the 
5 • 10 8 dark matter particles. 

The MGS consist of three simulation runs with different im- 
plementations of heating schemes: 

GO run: This simulation does not include any additional non- 
gravitational heating sources. As a consequence, it mostly re- 
produces the self-similar model expectations but drastically 
fails to reproduce the observed local galaxy cluster X-ray 
scaling relations, hence is not considered in the comparison 
analysis. 

PC run: This run employs a simplistic preheating model, raising 
the entropy of the gas particles to 200 keV cm 2 at z =4. In ad- 
dition to preheating, radiative cooling according to the cool- 
ing function of Sutherland & Dopita (1993) is implemented. 

FO run: The third simulation run includes no radiative cooling 
but incorporates a model of the energy input by supernova 
and AGN feedback computed by means of a semi-analytic 
model of galaxy formation, i.e. a gradual injection of energy 
into the surrounding ICM gas. 

The PC and FO runs represent two opposing heating 
schemes, one assuming an early preheating of the gas be- 
fore its accretion onto the cluster, the other incorporating non- 
gravitational heating as a relatively recent and ongoing process. 
Both runs are able to reproduce the observed local scaling re- 
lations, i.e. they predict the correct properties of the ICM for 
low-redshift clusters. However, owing to the contrasting heat- 
ing schemes, the two models make opposing predictions about 
the thermal evolution of the ICM with redshift. Table 4 gives an 
overview of the MGS evolution results fitted in the redshift range 
0<z<1.5. 

Table 4. Results on the redshift evolution of X-ray scaling rela- 
tions from the Millennium Gas Simulations (Short et al. 2010). 



GO run 


redshift evolution 


M-T 
Lx-T 
M-L x 


oc E(z)- 1 ■ (I + zf m±UMi 

oc E(z) 1 ■ (1 + z) 031±0M 
°c E(z)- 114 ■ (1 + z ) 0160±0 - 036 


PC run 


redshift evolution 


M-T 
Lx-T 
M-L x 


oc E(z)- 1 • (1 + z y u - m±UMl 
oc E(z) 1 ■ (I + z)- [ J1±0A6 
oc E(zT V4 ■ (1 + z)°- 594±om3 


FO run 


redshift evolution 


M-T 
Lx-T 
M-L x 


ocE(zr l .(i + z f«™ 

oc E(Z) 1 .(1 +Z )0 ^0.05 
oc E(z)- 114 ■ (1 + Z )-0.494±0.010 



Depending on the selection criteria of the cluster sample, the 
adopted fitting scheme, and since the fitted slope and normaliza- 
tion of scaling relations are not independent of each other, the 
measured normalization of local scaling relations differs by up 
to a factor of two (see e.g. the results of Arnaud et al. (2005) 




0.2 0.4 0.6 0.8 1 -L2 ti 

z 

Fig. 6. Redshift evolution of the M-T relation. Continuous red 
line and light grey confidence area: observed evolution. Green- 
dashed line and dark grey confidence area: MGS FO run. Blue- 
dashed line and dark grey confidence area: MGS PC run. Black- 
dashed line: self-similar prediction. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

z 

Fig. 7. Redshift evolution of the Lx-T relation. Red-dashed line: 
Estimated mean bias of the combined sample rescaled to remove 
the effects of bias in the local scaling relation. Additional lines 
and confidence areas have the same meaning as in Fig. 6. 

and Pratt et al. (2009)). Therefore, instead of directly comparing 
the redshift-dependent normalization derived from the MGS to 
the results of this work, we do not take into account the differ- 
ent local normalizations but only the evolution with respect to 
the local values. Fig. 6, 7, and 8 show the redshift evolution of 
the normalization of the M-T, L-T, and M-L^ relations with re- 
spect to its local value for both the MGS PC and FO runs and 
the observational results deduced in this work. We note that in 
contrast to this work, Short et al. (2010) assumed the self-similar 
evolution model and then fitted the observed difference from this 
model as powers of (1 + z). 

It is clearly visible that the FO simulation provides no good 
description of the observed evolution for all three relations. In 
contrast, the PC run shows good agreement with the observa- 
tions. For the M-T and M-Ly relations, the PC results are fully 
consistent with ours within the errors. The prediction for the Lx- 
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Fig. 8. Redshift evolution of the M-L^ relation. Lines and con- 
fidence areas have the same meaning as in Fig. 7. 



T relation shows slight deviations from our result at the 2<x level 
for z < 1.2. This difference can partly be attributed to the differ- 
ent functional form assumed for the evolution fit in Short et al. 
(2010). 

Care has to be taken when interpreting the constraints on the 
ICM thermal history deduced from this finding because none of 
the MGS runs provide a complete model of the necessary ICM 
heating and cooling processes. As an example, the PC run does 
not include ongoing heating, which is known to be of crucial im- 
portance to balance the cooling in cluster cores. Such an addition 
of mild ongoing feedback to the model is likely to bring the pre- 
dicted evolution for the L^-T relation in even closer agreement 
with the observed trend (see Fig. 7). Nevertheless, the observed 
evolution of X-ray scaling relations provides strong evidence of 
the preheating scenario. A further refinement of the simulations, 
e.g. the combination of preheating with mild ongoing feedback, 
is desirable to permit more detailed comparisons between the 
observed and the simulated evolutions. 

5.2. Implications for the eROSITA cluster survey 

eROSITA is the main instrument on the Spektrum-Rontgen- 
Gamma mission scheduled for launch in 2012 and will carry out 
a new all-sky X-ray survey in the energy range from 0.1 to 10 
keV (see Predehl et al. (2010) and Predehl et al. (2007)). One of 
the main science goals of the project is to provide a sample of 
~ 100000 X-ray selected galaxy clusters. A cluster catalogue of 
this size is necessary to test cosmological models to higher accu- 
racy and place reliable constraints on cosmological parameters 
such as the dark energy equation of state (Haiman et al. 2005). 

The achievable accuracy of the parameter constraints de- 
rived from the eROSITA all-sky survey depends on the number 
of cluster detections and the knowledge of X-ray scaling rela- 
tions up to high redshift. Whether a cluster can be detected by 
eROSITA is mainly determined by the cluster's soft band X-ray 
luminosity and its distance. Since the mean cluster luminosity 
at a given redshift depends on the evolution of X-ray scaling re- 
lations, the results of this work have a direct influence on the 
expected number of eROSITA clusters to be detected at high 
redshift and therefore on the expected cosmological constrain- 
ing power of the mission. 



The expected redshift distribution and total number of 
eROSITA clusters was recalculated by means of the cluster 
counter presented in chapter 8 of Miihlegger (2010). The cluster 
counter provides an estimate of the number of eROSITA clus- 
ter detections based on various simplifying assumptions. The 
assumed criterion for a cluster detection is whether the num- 
ber of photons detected from the observed system exceeds the 
count limit cn m . Until now, the eROSITA count limit has not been 
known in detail owing to its dependence on as of now unspeci- 
fied instrumental parameters such as the eROSITA point spread 
function, and furthermore its significant dependence on the to- 
tal X-ray background, which varies with sky position, exposure 
time and other conditions. A constant detection count limit of 
ciim = 100 is currently used as a conservative estimate. 

The number of clusters in each redshift and luminosity bin 
was set according to the luminosity function used in Sect. 2.4. 
This luminosity function is based on a cluster mass function re- 
sulting from a cosmological model consistent with the seven- 
year WMAP results (Komatsu et al. 2010) and was converted 
into a luminosity function by means of the L-M relation and the 
evolution results of Sect. 3.2. Further necessary input data for 
the cluster counter are an all-sky exposure map for the special 
geometry of the eROSITA survey and an all-sky map of Galactic 
neutral hydrogen. The ICM temperature was set according to the 
(ROSAT band) Lx-T relation of Pratt09 and the observed evolu- 
tion of Sect. 3.2. 

Thereafter, the eROSITA count rate at each position in the 
sky, cluster luminosity, and redshift is calculated by means of 
XSPEC, assuming an absorbed Mekal model with the param- 
eters z, Lx, T, and tin- The count rate is then converted into a 
number of detected photons by multiplying it with the exposure 
value at the respective sky position. If the number of detected 
photons for a given set of parameters exceeds ci; m , the number of 
clusters determined by the luminosity function for the given red- 
shift and luminosity is added to the cluster number in the current 
redshift bin and sky position. The resulting redshift distribution 
and total number of detected clusters is presented in the follow- 
ing section for three L^-T evolution models: The self-similar 
model, i.e. positive evolution, the no-evolution scenario assumed 
in Miihlegger (2010), and the slightly negative evolution found 
in this work. 

Table 5 shows the total number of expected cluster detec- 
tions with eROSITA under the assumption of a count limit 
ciim = 100. In addition to the all-sky expectations, the number 
of "extragalactic" detections refers to clusters with Galactic lat- 
itude |£>| >20°. The total number of achievable cluster detections 
is expected to be closer to this last number, since the high col- 
umn density of the absorbing Galactic interstellar medium and 
high density of other X-ray and stellar sources in the Galactic 
plane make cluster detections in this area challenging. For the 
extragalactic area, the number of expected clusters with z > 
0.8, 1, 1.2, 1.4, and 1.6 is also listed. 

As can be seen in Fig. 9, which shows the number of cluster de- 
tections with redshift within the redshift interval [z - 0.005; z + 
0.005] and Fig. 10, the achievable number of cluster detections 
with redshift > z, a change in the assumed evolution model has 
a significant influence on the number of eROSITA cluster de- 
tections. However, the total number of clusters that should be 
detected is hardly affected because the majority of the eROSITA 
clusters are expected to be discovered at low redshifts, where 
the scaling relation evolution is of little importance. As an ex- 
ample, the total number of clusters for the extragalactic area 
is decreased by ~ 5% when changing from the previously as- 
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Table 5. Number of expected cluster detections with eROSITA 
assuming self-similar, no, or slightly negative evolution of the 
Lx-T relation. The rows labeled "extragalactic" refer to a 
Galactic latitude \b\ > 20°. A count limit of ci; m = 100 was as- 
sumed. 
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N :>1 




N z >1.4 




all sky 
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extragal. 


97195 
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Fig. 9. Achievable number of cluster detections with redshift 
within the redshift interval [z - 0.005; z + 0.005] for eROSITA. 
Continuous lines: All sky. Dashed lines: Extragalactic (\b\>20°). 
Black: Self-similar Lx-T evolution (oc E(z) +l ). Green: No Lx-T 
evolution. Red: Best fit Lr-T evolution 



sumed no-evolution scenario to the negative evolution found in 
this work. 

The situation is clearly different for distant cluster detec- 
tions. Owing to the smaller number of luminous high-redshift 
clusters in the negative evolution scenario, the number of ex- 
pected detections is smaller by ~ 18% at z = 0.8 and by ~28% at 
Z= 1.2 relative to the no-evolution scenario. However, even with 
these corrected cluster expectations, eROSITA is still likely to 
increase the sample of known z > 0.8-clusters by about a factor 
of 50 from present numbers. 

6. Summary and conclusions 

The main goal of this study has been to to investigate the redshift 
evolution of galaxy cluster X-ray scaling relations by means of 
a combined cluster sample. To this end, a cluster sample was 
compiled from both recent publications and newly discovered 
clusters provided by the XMM-Newton Distant Cluster Project 
(XDCP). Our gathered sample of recently discovered distant 
clusters has allowed tighter constraints to be made on scaling 
relation evolution than possible in previous studies. 
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Fig. 10. Achievable number of cluster detections with redshift 
>z for eROSITA. The lines have the same meaning as in Fig. 9. 



The definition of cluster observables slightly differs between 
the individual publications used to compile the combined cluster 
sample. A homogenization scheme that accounts for these dif- 
ferent analysis schemes was therefore applied to the subsamples. 
In detail, the definition of cluster radii, either relying on a fixed 
density contrast A at the cluster redshift or a redshift-dependent 
average overdensity A ; and the different values for the density 
contrast were corrected for. Furthermore, the applied corrections 
take into account the slight differences in the assumed cosmolog- 
ical parameters, in the energy band for which Lx is given. 

On the basis of data for clusters at z < 0.3 included in the 
combined sample, local M-T, Lx-T, and M-L -relations were 
fitted. The derived relations generally agree well with the pub- 
lished results of Pratt09 and show deviations from the self- 
similar model similar to those found in earlier studies, e.g. a 
steeper Lr-T relation. 

Typical X-ray selected cluster samples can be approximated 
to be flux-limited. Various selection biases complicate the anal- 
ysis of these flux-limited samples and are important to the in- 
terpretation of results on scaling relations and their evolution. 
The bias inherent to the sample used in this work was estimated 
by means of comparable cluster samples selected from a simu- 
lated cluster population. In detail, selection biases were found to 
raise the measured normalization of the local Lx-T relation and 
decrease the apparent scatter about the mean relation for high-z 
clusters. For the Lx-T relation, the bias appears to generate a 
positive evolution, whereas for the M-Lx relation it generates 
the opposite. 

The redshift evolution of X-ray scaling relations was investi- 
gated in the redshift range <z < 1.46. Throughout this redshift 
range, no significant variation in the slope of the relations was 
found. The normalization, however, evolves with redshift for all 
three examined relations. For the M-T relation, the measured 
evolution is ocis^)" 104 * 07 , consistent with the self-similar pre- 
diction. The results for the L x -T (oc E(zy 02r -»£) and M-Lx 
relation (oc E(zy 093 -°"), however, differ significantly from the 
model predictions. As for the local relations, these deviations in- 
dicate that the influence of non-gravitational ICM heating and 
cooling is not negligible. The inconsistent results of recent stud- 
ies have been found to be caused by limited sample sizes and red- 
shift ranges in combination with selection biases and to a lesser 
degree by the use of different local relations. 
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On the basis of our results and the local scaling of Pratt09 
and assuming that h = 0.70, the M-T, L-T, and M-L x relations 
for cluster properties within rjoo have the explicit form 



M = (0.291+0.031)-(T[keV]) 



1.62±0.08 



£ (Z )-1-04±0.07 10 14 M 



(23) 



L x = (0.079+0.008)<r[keV]) 2 - 70±0 - 24 -£(z)" ' 23! »^10 44 erg s _1 , (24) 



M=(1.39±0.07)-(L x [10 44 erg s" 1 ]) 1 



-l n \0.54±0.03 



£(z)-°- 93 -o-10 14 M R -(25) 



This work once again highlights the importance of selection bi- 
ases to scaling relation studies. For distant cluster mass estimates 
based on the bolometric X-ray luminosity Lx, we recommend 
using Equ. 26, which is based on the bias-corrected local relation 
of Pratt09 and includes a correction for the estimated selection 
bias on the observed evolution, given by 



M = (1.64±0.07)-(L x [10 44 erg s" 1 ])' 



-1 -,-,0.52*0.03 



£(z)- a90 -o-10 14 M (26) 



To provide tighter constraints on scaling-relation evolution and 
improve the mass estimate of Equ. 26 in the future, a more homo- 
geneous, extensive distant cluster sample with a precisely known 
selection function is necessary. Such a sample would allow a 
more precise bias estimate and correction to be made than pos- 
sible for the sample used in this study. 

Comparing the observed evolution with predictions made by 
the Millennium Gas Simulations has allowed us to discriminate 
between different proposed scenarios and attempt a physical in- 
terpretation of the thermodynamic history of the ICM. The com- 
parison analysis strongly suggests an early preheating, i.e. an en- 
tropy increase for the gas particles before the infall of the ICM 
gas into the cluster potential well. 

The expected number of cluster detections for the upcom- 
ing eROSITA survey was recalculated taking into account the 
results of this work. In general, the total number of achievable 
detections is slightly lower than assumed before, while the num- 
ber of high redshift clusters to be detected shows a significant 
decrease. 

Future more detailed studies of the redshift evolution of X- 
ray scaling relations will be important to more tightly constrain 
the early thermodynamic history of the ICM and provide cali- 
brated mass-observable relations for upcoming large cluster sur- 
veys and their cosmological applications. 
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Table A.l. Clusters included in the combined sample and their X-ray properties within rsoo- T indicates the global ICM temperature, 
Lx the bolometric X-ray luminosity and M the (hydrostatic) cluster mass. Cluster names printed in boldface designate the additional 
z> 0.8 -clusters added to the sample (see Sect. 2.1). 



Cluster 


z 


Kl [K&\\ 


t r 1 r\44 ...» „ — 1 n 

L x [10 erg s J 




Source publication 


XMMXCSJ2215.9-1738 


1.46 




z - z - , -0.35 




Hilton et al. (2010) 


ISCSJ1438.1+3414 


1.41 




z - z '- , -0.53 




Brodwin et al. (2010) 


XMMUJ2235.3-2557 


1.39 


8.6![| 


10.0 + 0.8 


4.4 + 1 


Rosati et al. (2009) 


RX J0848+4453 


1.270 


2.9 + 0.8 


1.0 + 0.7 


1.2 ±0.9 


Ettori04 


RX J0849+4452 


1.260 


5.2 ± 1.6 


2.8 ± 0.2 


2.7 ± 1.4 


Ettori04 


RX J 1252-2927 


1.235 


5.2 + 0.7 


5.8 ± 1.0 


1.5 ±0.3 


Ettori04 


XLSSJ022303.0-043622 


1.22 


3.8!~ 


1.1 ±0.7 




Bremer et al. (2006) 


RXJ1053.7+5735 


1.134 


3.9 + 6.2 


2.2 + 0.7 




Hashimoto et al. (2004) 


SPT-CLJ2106-5844 


1.132 


8.5^-6 


73.1+5.3 




Foley et al. (2011) 


RX J09 10+5422 


1.106 


6.6 ± i.7 


2.8+0.3 


4.6 ±2.8 


Ettori04 


XMMUJ100750.5 +125818 


1.082 


5 7 +co 


1 ^ + °° 

i - J -0.5 




Schwope et al. (2010) 


SPTJ0546-5345 


1.0665 


7.5 + 1.7 


26.9 ± 1.7 




AnderssonlO 


3C186 


1.067 


- - r,_i_n 98 

j.jo_ 027 


13.5 + 0.7 


~ , i A 1 
'J J 1 TU.Al 

-P44 


Siemiginowska et al. (2010) 


XLSSJ022404.1-041330 


1.050 


4.1 + 0.9 


-) A 4-0 3 

J -0.2 


1 O 4-0 Q 

i - 3 -0.3 


Maughan et al. (2008) 


XLSSJ022709. 2-041800 


1.05 


3.7 ± 1.5 


1.9 + 0.2 


1.3 ±0.3 


Pacaud07 


C1J1415. 1+3612 


1.03 


5.7 ± 1.2 


10.1 ± 0.6 


3.0 ±0.9 


Maughan06 


SPTJ2341-5119 


0.9983 


8.0 ± 1.9 


24.5 ± 1.3 




AnderssonlO 


2XMMJ083026.2+524133 


0.99 


8.2 ± 0.9 


18.5 ± 3.6 


5.6 ± 1 


Lamer et al. (2008) 


XMMUJ1229.5+0151 


0.975 


6 4 +0 - 7 


8.8 + 1.5 




Santos et al. (2009) 


XMMUJ1230.3+1339 


0.975 


5 3 +0 - 7 

J - -0.6 


6.5+0.7 


3.1 ± 1.8 


Fassbender et al. (201 1) 


CU1429.0+4241a 


0.92 


6.2+1.5 


9.3+0.9 


3.5 ± 1.3 


Maughan06 


RCS23 19+0030 


0.904 


6 ± 1 


7.9 + 0.7 


2.8 ±0.3 


Hicks08 


CLJ1604+4304 


0.90 


2 51 +105 

-0.69 


2 + 0.4 




Lubin et al. (2004) 


RCS2320+0033 


0.901 


6.0 + 2 


5.9 + 0.5 


2.5 ±0.5 


Hicks08 


RCS23 19+0038 


0.900 


5.9 + 0.8 


16.2 + 0.6 


3.6 ±0.3 


Hicks08 


WGA1226+3333 


0.890 


11.2 + 2.2 


54.2 ± 0.8 


9.1 ±2.4 


Ettori04 


SPTJ0533-5005a 


0.8810 


4.0+1.9 


3.6 + 0.9 




AnderssonlO 


C1J1008.7+5342 


0.87 


3.6 + 0.8 


3.6 + 0.3 


1.5 ±0.6 


Maughan06 


RCS 1620+2929 


0.870 


3.9 + 1.2 


3.3 + 0.5 


2.1 ±0.4 


Hicks08 


RXJ1257.2+4738 


0.866 


3 6 +2 ' 9 


2 +29 




Ulmer et al. (2009) 


ClJ1559.1+6353a 


0.85 


4.1 ± 1.4 


2.5+0.3 


1.6+1.1 


Maughan06 


XLSSJ022738.3-031758 


0.84 


3.3 + 1.1 


4.1 ±0.4 


0.9 ± 0.2 


Pacaud07 


RX J0152-1357N 


0.835 


6.0+1.1 


10.2 + 0.6 


2.4 ± 0.7 


Ettori04 


RXJ0152-1357S 


0.830 


6.9 + 2.9 


7.5 ±0.4 


3.2 ± 1.4 


Ettori04 


MS 1054-0321 


0.830 


10.2+1.0 


28.4 ±3.0 


19.1+3.5 


Ettori04 


C1GJ1056-0337 


0.826 


9.2+1.5 


49.1 ± 1.2 




OHara07 


RXJ1821.6+6827 


0.816 


4 7 +1 ' 2 

-0 7 


io.4 + ;-3 




Gioia et al. (2004) 


RX J1716+6708 


0.813 


6.8+ i.O 


13.6 + i.O 


4.0 + 0.8 


Ettori04 


RX J1350+6007 


0.810 


4.6 ± 0.7 


4.2 + 0.4 


1.4 + 0.4 


Ettori04 


RX J1317+2911 


0.805 


4.1 ± 1.2 


1.1 ±0.1 


1.4 + 0.6 


Ettori04 


XLSSJ022210.7-024048 


0.79 


3.9 + 2.8 


1.4 ±0.2 


1.8 + 0.4 


Pacaud07 


ClJ1103.6+3555a 


0.78 


6.0 + 0.9 


4.5 ±0.2 


2.9 + 0.7 


Maughan06 


RCS2318+0034 


0.78 


5.8 + 1.2 


7.7 ±0.8 


12.9 + 2.0 


Hicks08 


MSI 137+6625 


0.782 


6.9 + 0.5 


15.2 ±0.4 


4.7 ± 0.6 


Ettori04 


SPTJ2337-5942 


0.7814 


8.9 + 2.0 


41.3 ±2.3 




AnderssonlO 


RCS0224-0002 


0.778 


5.1 ± 1.2 


4.1 ±0.5 


5.0 + 0.7 


Hicks08 


XLSSJ022532.2-035511 


0.77 


2.8 + 0.8 


2.1 ±0.2 


0.9 ± 0.2 


Pacaud07 


SPTJ0528-5300a 


0.7648 


5.2 + 3.5 


6.6 ±0.7 




AnderssonlO 


RCS 1107-0523 


0.735 


4.2 ± 0.6 


3.3 ±0.3 


1.8 ±0.2 


Hicks08 


RX J2302+0844 


0.734 


6.6+1.5 


5.3 ±0.2 


3.2 ±0.7 


Ettori04 


RXJ1 113-2615 


0.730 


5.6 + 0.8 


4.4 ±0.8 


3.2 ±0.7 


Ettori04 


C1J1113. 1-2615 


0.72 


4.7 + 0.9 


3.7 ±0.3 


2.6 ±0.8 


Maughan06 
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Cluster 


7 


kT[keV] 


LxtlO^ergs- 1 ] 


M[10 14 MJ 




Source publication 


ClJ1342.9+2828a 


0.71 


3.7 + 0.5 


3.3 ±0.1 


1.8 ±0.4 


Maughan06 


WJ 1342.8+4028 


0.70 


3.5+0.3 


3.8 ±0.8 


1.3 ±0.2 


Kotov05 


RX J1221+4918 


0.700 


7.5 + 0.7 


12.7 ±0.4 


5.5 ±0.9 


Ettori04 


C1GJ0744+3927 


0.686 


9.6 ± 0.9 


46.7 ± 0.7 




OHara07 


C1GJ1419+5326 


0.640 


4.1 ±0.8 


6.9 ±0.2 




OHara07 


RX J0542-4100 


0.634 


7.9+1.0 


11.6+1.3 


3.9 ±0.8 


Ettori04 


RCS 1419+5326 


0.62 


4.6 ± 0.4 


8.4 ± 0.5 


2.6 ±0.147 


Hicks08 


RXJ1334.3+5030 


0.62 


4.6 ± 0.4 


7.7 ± 1.5 


2.8 ±0.5 


Kotov05 


C1J0046.3+8530 


0.62 


4.4 ± 0.5 


3.8 ±0.2 


2.1 ±0.4 


Maughan06 


XLSSJ022457. 1-034856 


0.61 


3.2 ± 0.4 


3.6 ±0.2 


1.2 ±0.3 


Pacaud07 


SPTJ0559-5249 


0.6112 


7.7 ± 1.1 


14.1 ±0.9 




AnderssonlO 


RXJ1120.1+4318 


0.60 


4.9 ± 0.3 


13.4 ±2.7 


4.7 ± 1.2 


Kotov05 


C1GJ0647+7015 


0.584 


15.0 + 3.8 


49.9 ± 1.1 




OHara07 


MS2053-0449 


0.583 


5.5 + 0.5 


5.3 ± 1.0 


3.1 ±0.5 


Ettori04 


SPTJ233 1-5051 


0.5707 


5.9 + 1.3 


16.2 ±0.7 




AnderssonlO 


RX J0848+4456 


0.570 


3.2 + 0.3 


1.2 ±0.4 


1.4 ±0.3 


Ettori04 


C1GJ2 129-0741 


0.570 


11.8 + 2.8 


39.4 ± 1.3 




OHara07 


RXJ1 121+2326 


0.562 


4.6 + 0.5 


5.4 ± 0.2 


5.1 ± 1.8 


Ettori04 


C1GJ0717+3745 


0.548 


11.5 + 0.7 


121.8 ±2.2 




OHara07 


C1GJ 1354-0221 


0.546 


4.1 ±0.8 


5.9 ± 0.3 




OHara07 


C1GJ1423+2404 


0.545 


5.4 + 0.2 


37.2 ±0.2 




OHara07 


C1GJ 1149+2223 


0.544 


9.8 ±0.8 


64.1 ± 1.5 




OHara07 


MS0016+1609 


0.541 


10.0 ±0.5 


52.0 ±7.2 


8.8 ±0.7 


Ettori04 


MS045 1-0305 


0.539 


8.0 ±0.3 


50.2 ±7.9 


7.1 ±0.5 


Ettori04 


RX J 1525 +0957 


0.516 


5.1 ±0.5 


6.6 ± 0.2 


2.8 ±0.5 


Ettori04 


RX J0505.3+2849 


0.51 


2.5 ±0.4 


1.6 ±0.3 


1.4 ±0.2 


Kotov05 


XLSSJ022206.7-030314 


0.49 


3.6 ±0.6 


3.1 ±0.2 


1.6 ±0.3 


Pacaud07 


XLSSJ022357.4-043517 


0.49 


2.2 ±0.9 


0.46 ± 0.04 


0.5 ±0.1 


Pacaud07 


CL 1641+4001 


0.464 


5.1 ±0.8 


3.0 ± 0.2 




Branchesi07 


SPTJ0509-5342 


0.4626 


7.0 ± 1.4 


13.0 ±0.4 




AnderssonlO 


C1GJ1621+3810 


0.461 


6.8 ±0.6 


20.2 ± 0.2 




OHara07 


3C 295 


0.460 


4.3 ±0.3 


13.7 ±0.2 


2.0 ±0.2 


Ettori04 


RX J1701+6421 


0.453 


4.5 ±0.4 


5.8 ± 0.5 


1.3 ±0.1 


Ettori04 


RXJ1347-1145 


0.451 


10.3 ± 0.6 


114.8 ±0.6 


7.8 ±0.7 


Ettori04 


C1GJ0329-0212 


0.450 


5.9 ±0.2 


30.5 ± 0.2 




OHara07 


MACSJ0417.5-1154 


0.440 


9.4 ±0.7 


128.2 ±2.0 




OHara07 


RXCJ1206.2-0848 


0.440 


11.4 ±0.9 


53.0 ±0.6 




OHara07 


XLSSJ022145.2-034617 


0.43 


4.8 ±0.6 


6.6 ± 0.2 


2.8 ±0.6 


Pacaud07 


MS 162 1+2640 


0.426 


6.8 ±0.9 


10.3 ±0.3 


4.0 ±0.8 


Ettori04 


MS 0302.5+1717 


0.42 


4.5 ±0.5 


4.4 ± 0.9 


2.2 ±0.6 


Kotov05 


MS0302+1658 


0.424 


3.8 ±0.9 


6.3 ± 0.2 


2.1 ±0.5 


Ettori04 


SPTJ0551 -5709b 


0.4230 


4.1 ±0.9 


5.8 ± 0.6 




AnderssonlO 


RXCJ2228.6+2036 


0.412 


8.1 ±0.5 


35.2 ±0.7 




OHara07 


RX J1416+4446 


0.400 


3.7 ±0.2 


5.0 ± 0.2 


1.2 ±0.1 


Ettori04 


Zw CI 0024.0+1652 


0.394 


4.4 ±0.5 


4.2 ± 0.2 




Branchesi07 


RXCJ0949.8+1707 


0.383 


7.8 ±0.7 


31.3 ±0.5 




OHara07 


Zw CI 1953 


0.374 


7.6 ±0.5 


26.2 ± 0.4 




OHara07 


Abell 370 


0.373 


8.7 ±0.5 


22.4 ± 0.4 




OHara07 


RXJ1532.9+3021 


0.362 


6.1 ±0.3 


43.8 ±0.2 




OHara07 


RXCJ0404.6+1109 


0.355 


5.6 ±0.8 


9.8 ± 0.6 




OHara07 


XLSSJ022722.4-032144 


0.33 


2.4 ±0.5 


0.7 ±0.1 


0.9 ±0.2 


Pacaud07 


Abell 1722 


0.328 


9.1 ± 1.5 


13.7 ±0.3 




OHara07 


ZwCl 1358+6245 


0.327 


9.1 ±0.9 


19.4 ±0.3 




OHara07 


XLSSJ022402.0-050525 


0.32 


2.0 ±0.7 


0.14 ±0.02 


0.7 ±0.1 


Pacaud07 


Abell 1995 


0.318 


8.1 ± 1.0 


17.1 ±0.2 




OHara07 


MS2137. 3-2353 


0.313 


5.0 ± 0.2 


32.1 ±0.2 




OHara07 
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Cluster 


7 


kT[keW] 


LxtlO^ergs- 1 ] 


M[10 14 MJ 




Source publication 


Abell 1300 


0.3075 


9.2 + 0.4 


18.0 ± 1.5 


5.2 + 3.0 


Zhang07 


Abell 2744 


0.3066 


10.1+0.3 


22.12 ± 1.7 


7.4 + 2.9 


Zhang07 


RXCJ2245. 0+2637 


0.304 


5.9 + 0.3 


18.3 ±0.1 




OHara07 


MS 1008.1-1224 


0.302 


6.0 + 0.4 


11.3 ±0.4 




Branchesi07 


Abell 781 


0.298 


6.5 + 0.5 


6.3 ±1 


4.5 ± 1.3 


Zhang08 


RXC J2308.3-0211 


0.297 


7.6 + 0.7 


12.1 ± 1.3 


7.4 + 2.3 


Zhang08 


Abell 2537 


0.2966 


7.6 + 0.9 


15.8 ± 1.0 


7.2+1.1 


Mantz09 


1ES 0657-558 


0.2965 


11.7 + 0.2 


65.2 ±0.9 


22.8 ± 2.8 


Mantz09 


RXC J0658.5-5556 


0.296 


10.7 ± 0.4 


49.5 ± 2.4 


11.0 ± 5.3 


Zhang08 


RXC J05 16.7-5430 


0.294 


6.7+0.5 


10.9 ± 1.5 


6.4 ± 1.8 


Zhang08 


XLSSJ022803.4-045103 


0.29 


2.8 + 0.6 


0.50 ± 0.04 


1.4 ±0.1 


Pacaud07 


RXC J0043 .4-2037 


0.292 


7.0 + 0.4 


10.5 ± 1 


4.8 ± 1.4 


Zhang08 


Zw CI 3146 


0.2906 


8.4 + 0.4 


49.1 ± 1.8 


9.4 ± 1.2 


Mantz09 


Abell 611 


0.288 


8.9 + 0.7 


18.9 + 0.3 




OHara07 


RXC J0528.9-3927 


0.284 


6.6 + 0.5 


14.5 ± 1.5 


6.4 ± 1.9 


Zhang08 


RXC J0232.2-4420 


0.284 


6.6 + 0.3 


18.5 ± 1.4 


8.4 ±2.5 


Zhang08 


RXCJ0437. 1+0043 


0.2842 


5.1+0.3 


6.2 ± 0.7 


6.1 ±2.2 


Zhang07 


Abell 697 


0.282 


10.9+ 1.1 


41.9 ±2.3 


17.1 ±2.9 


Mantz09 


Abell 1758 


0.280 


7.9 + 0.2 


10.5 ±0.7 


11.2 ±3.4 


Zhang08 


RXJ201 1.3-5725 


0.2786 


3.2 + 0.3 


6.0 ± 0.3 


3.3 ±0.7 


Mantz09 


RXC J0532.9-3701 


0.275 


7.7+0.6 


12.7 + 1.2 


5.4 ± 1.6 


Zhang08 


RXC J2337.6+0016 


0.275 


7.5+0.5 


10.0 + 0.9 


11.0 ± 3.2 


Zhang08 


RXCJ0303.7-7752 


0.2742 


8.2 + 0.5 


12.9+1.3 


7.7 ±2.3 


Zhang07 


XLSSJ022524.7-044039 


0.26 


2.0 + 0.2 


0.49 ± 0.02 


0.6 ±0.1 


Pacaud07 


MS 1006.01202 


0.261 


6.1 ±0.4 


10.8 + 0.2 




OHara07 


RXC J0307. 0-2840 


0.258 


7.1 ±0.4 


13.1 ± 1.2 


5.5 ±2.0 


Zhang08 


Zw CI 7160 


0.258 


4.7+0.1 


17.0 ±0.8 


2.4 ±0.7 


Zhang08 


MS J1455.0+2232 


0.2578 


4.5 ±0.2 


20.7 ± 0.6 


6.2 ± 1.0 


Mantz09 


Abell 68 


0.255 


7.3 ±0.3 


11.4+1 


6.5 ± 1.9 


Zhang08 


Abell 1835 


0.253 


8.4 ±0.3 


53.2 ± 1.7 


8.0 ±2.3 


Zhang08 


Abell 521 


0.2475 


6.2 ±0.3 


15.7 ±0.7 


11.4+ 1.7 


Mantz09 


Abell 2125 


0.246 


3.4 ±0.2 


1.9 ±0.1 




Branchesi07 


RX J0439.0+0715 


0.2443 


6.6 ±0.5 


18.4 ±0.7 


7.4 ± 1.0 


Mantz09 


RXC J2129.6+0005 


0.235 


6.3 ±0.2 


14.3 ±0.9 


4.3 ± 1.3 


Zhang08 


Zw CI 2089 


0.2347 


6.6 + 1.5 


13.2 ±0.7 


3.1 ±0.4 


Mantz09 


Abell 2390 


0.233 


11.6 ±0.6 


40.9 ± 2.7 


7.7 ±2.3 


Zhang08 


Abell 2667 


0.230 


7.0 ±0.3 


21.6 ± 1.3 


6.0 ± 1.7 


Zhang08 


Abell 267 


0.230 


6.2 ±0.4 


7.7 ± 0.7 


4.3 ± 1.3 


Zhang08 


Abell 21 11 


0.229 


6.5 ±0.7 


11.3 ± 0.7 


7.8 ± 1.9 


Mantz09 


Zw CI 5247 


0.229 


5.3 ± 1.1 


8.8 ± 0.6 


8.2 ± 1.8 


Mantz09 


Abell 1763 


0.228 


5.8 + 0.3 


16.5 ± 1.5 


5.0 ± 1.5 


Zhang08 


RX J0220.9-3829 


0.228 


5.2 + 0.5 


8.0 ± 0.4 


4.4+1.1 


Mantz09 


Abell 2219 


0.2281 


10.9 + 0.5 


45.1+2.3 


18.9 ±2.5 


Mantz09 


Abell 1682 


0.226 


7.0 + 2.1 


15.3 + 1.7 


12.4 ±3.2 


Mantz09 


RX J0638.7-5358 


0.2266 


9.5 ± 1.0 


30.6+1.6 


10.3 ± 1.4 


Mantz09 


Abell 2261 


0.224 


6.6 + 0.6 


13.9 + 2.2 


6.0 ± 1.7 


Zhang08 


RX J0237.4-2630 


0.2216 


6.7 + 1.3 


12.8 + 0.7 


5.6 + 1.1 


Mantz09 


RXJ0304. 1-3656 


0.2192 


6.3 ±0.8 


6.7 ± 0.4 


4.4 ±0.9 


Mantz09 


Abell 773 


0.217 


8.3 ±0.4 


20.9+1.6 


8.3 ±2.5 


Zhang08 


RX J 1504. 1-0248 


0.2153 


8.0 ±0.4 


69.4 ± 2.5 


11.0 ± 1.4 


Mantz09 


Zw CI 2701 


0.214 


6.8 ±0.5 


10.4 ±0.5 


4.0 ±0.7 


Mantz09 


Abell 1423 


0.213 


5.8 ±0.6 


13.2 ±0.9 


8.7 ±2.0 


Mantz09 


Abell 209 


0.209 


7.1 ±0.3 


13.3 ± 1.1 


5.3 ± 1.7 


Zhang08 


RX J0439.0+0520 


0.208 


5.0 ±0.5 


8.7 ± 0.2 


2.7 ±0.5 


Mantz09 


Abell 963 


0.206 


6.5 ±0.2 


11.4 + 0.9 


5.2 ± 1.5 


Zhang08 


Abell 520 


0.203 


7.2 ± 0.2 


20.1+0.7 


11.9 ± 1.6 


Mantz09 
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Cluster 


z 


kT[keW] 


LxtlO^ergs- 1 ] 


M[10 I4 M H ] 


Source publication 


Abell2163 


0.203 


12.3 + 0.9 


88.1+3.4 


38.5+5.0 


Mantz09 


Abell 115 


0.197 


6.2 + 0.1 


14.3 + 1.1 


4.2+1.1 


Zhang08 


Abell 383 


0.187 


4.7 + 0.2 


8.1+0.5 


3.2 + 0.9 


Zhang08 


Abell 1689 


0.184 


8.5 + 0.2 


28.4 ± 1 


10.3 + 3.0 


Zhang08 


RXCJ131 1.4-0120 


0.1832 


8.9 + 0.1 


36.1+0.1 




Pratt09 


Abell 665 


0.1818 


8.0 + 0.2 


21.7 + 2.0 


12.7+ 1.8 


Mantz09 


Abell 2218 


0.176 


6.6 + 0.3 


11.1+0.8 


4.2+1.3 


Zhang08 


Abell 1914 


0.171 


8.8 + 0.3 


21.7 ± 1.1 


16.8 + 4.9 


Zhang08 


RXC J0958.3-1103 


0.167 


5.8 + 0.3 


9.1+0.9 


3.3 + 1.0 


Zhang08 


RXC J0645.4-5413 


0.164 


7.6 + 0.3 


17.8 ± 1.6 


6.6 ± 2.0 


Zhang08 


Abell 901 


0.163 


3.2 + 0.2 


16 + 0.5 


3.2+1.0 


Zhang08 


Abell 907 


0.1603 


5.96 ± 0.08 




4.7 ± 0.4 


Vikhlinin06 


RXC J2014.8-2430 


0.1538 


4.8 + 0.1 


21.1 ±0.1 




Pratt09 


RXC J0945.4-0839 


0.153 


5.3 + 0.5 


4.2 ± 0.6 


6.7 ± 2.0 


Zhang08 


Abell 2204 


0.152 


7.6 + 0.2 


33.9+1.2 


5.3 + 1.5 


Zhang08 


RXC J2234.5-3744 


0.151 


7.5 + 0.4 


10.9+1.1 


8.1 ±2.4 


Zhang08 


RXC J2217.7-3543 


0.1486 


4.9 + 0.1 


6.12 + 0.03 




Pratt09 


RXC J0547.6-3152 


0.148 


6.0 + 0.2 


7.2 + 0.6 


5.8 ± 1.7 


Zhang08 


RXC J2048. 1-1750 


0.1475 


4.7 + 0.1 


5.13 + 0.03 




Pratt09 


Abell 1413 


0.143 


6.6 + 0.1 


13.8 + 0.8 


5.4+1.6 


Zhang08 


XLSSJ022540.6-031121 


0.14 


3.5 + 0.6 


1.0 + 0.1 


2.4 + 0.1 


Pacaud07 


RXC J2218.6-3853 


0.141 


6.2 + 0.2 


6.6 ± 0.5 


4.8 ± 1.4 


Zhang08 


RXC J0020.7-2542 


0.1410 


5.7 + 0.1 


6.52 ± 0.04 




Pratt09 


RXC J0605. 8-3518 


0.1392 


4.6 + 0.1 


9.54 ± 0.04 




Pratt09 


Abell 1068 


0.1386 


3.9 + 0.1 


9.7 + 0.1 


3.7+0.6 


Mantz09 


RXC J1044.5-0704 


0.1342 


3.41+0.03 


7.42 ± 0.02 




Pratt09 


MS Jllll. 8 


0.1306 


5.8 + 0.2 


6.8 ± 0.6 


4.5 ±0.7 


Mantz09 


RXC J1516. 5-0056 


0.1198 


3.6 + 0.1 


2.31+0.02 




Pratt09 


RXC Jl 141.4-1216 


0.1195 


3.31+0.03 


3.75 + 0.01 




Pratt09 


RXC J2149. 1-3041 


0.1184 


3.26 ± 0.04 


3.56 ± 0.02 




Pratt09 


RXC J1516. 3+0005 


0.1181 


4.5 + 0.1 


4.12 + 0.02 




Pratt09 


RXC J0145.0-5300 


0.1168 


5.5 + 0.1 


5.00 ± 0.03 




Pratt09 


RXC J0616.8-4748 


0.1164 


4.2 + 0.1 


2.38 ± 0.02 




Pratt09 


RXC J0006.0-3443 


0.1147 


5.0 + 0.2 


4.1+0.1 




Pratt09 


Abell 2034 


0.113 


7.2 + 0.3 


9.5+1.0 


6.7 ± 1.0 


Mantz09 


RXC J0049.4-2931 


0.1084 


3.1+0.1 


1.78 + 0.02 




Pratt09 


PKS 0745-191 


0.1028 


8.0 + 0.3 




7.3+0.8 


Arnaud05 


RXC J021 1.4-4017 


0.1008 


2.1+0.1 


0.81+0.01 




Pratt09 


Abell 2244 


0.0989 


5.4 + 0.1 


10.7 ± 1.0 


6.2+1.1 


Mantz09 


RXC J2319.6-7313 


0.0984 


2.22 ± 0.03 


2.00 ± 0.02 




Pratt09 


Abell 3921 


0.094 


5.1+0.2 


6.2 + 0.6 


5.4 + 0.9 


Mantz09 


RXC J0003.8+0203 


0.0924 


3.9 + 0.1 


1.88 + 0.01 




Pratt09 


Abell 2142 


0.0904 


10.0 + 0.3 


34.8+1.7 


13.9 + 2.1 


Mantz09 


Abell 478 


0.0881 


8.0 + 0.3 


33.4 + 3.0 


10.1 ± 1.6 


Mantz09 


Abell 2597 


0.0852 


3.6 + 0.1 


7.2 + 0.7 


2.9 + 0.5 


Mantz09 


RXC J 1302. 8-0230 


0.0847 


3.0 + 0.1 


1.38 + 0.01 




Pratt09 


RXC J0821.8+0112 


0.0822 


2.7 + 0.1 


0.77 ± 0.01 




Pratt09 


Abell 2255 


0.0809 


6.4 + 0.2 


6.5 + 0.7 


5.9 ± 1.0 


Mantz09 


RXC J2129.8-5048 


0.0796 


3.8 + 0.2 


1.46 + 0.02 




Pratt09 


RXC J1236. 7-3354 


0.0796 


2.7 + 0.1 


1.03 + 0.01 




Pratt09 


Abell 2029 


0.0779 


8.2 + 0.2 


27.0 ± 2.6 


9.3 ± 1.4 


Mantz09 


Abell 3112 


0.0752 


4.3 + 0.1 


8.0 + 0.7 


4.1+0.6 


Mantz09 


Abell 401 


0.0743 


7.7 + 0.3 


16.8+1.0 


10.1 ± 1.6 


Mantz09 


Abell 1795 


0.0622 


6.1+0.1 


13.2 + 0.9 


5.5 + 0.8 


Mantz09 


RXC J0225. 1-2928 


0.0604 


2.5 + 0.2 


0.51+0.01 




Pratt09 


RXC J0345.7-4112 


0.0603 


2.19 ± 0.04 


0.77 ± 0.01 




Pratt09 
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Cluster 


z 


itrfkeV] 


LxtlO^ergs- 1 ] 


M[10 14 MJ 


Source publication 


Abell 3266 


0.0602 


8.6 ±0.2 


12.8+0.8 


9.2 ± 1.4 


Mantz09 


Abell 1991 


0.0586 


2.7 + 0.1 




1.2 ±0.1 


Arnaud05 


Abell 2256 


0.0581 


6.9 + 0.2 


10.7+0.9 


7.2 ± 1.0 


Mantz09 


RXC J2157.4-0747 


0.0579 


2.5 + 0.1 


0.45 + 0.01 




Pratt09 


Abell 133 


0.0569 


4.1+0.1 




3.3 ± 0.4 


Vikhlinin06 


RXC J2023.0-2056 


0.0564 


2.7 + 0.1 


0.61+0.01 




Pratt09 


Abell 85 


0.0557 


6.5 + 0.1 


12.8 ±0.9 


7.2 ± 1.0 


Mantz09 


Abell 3667 


0.0557 


6.3 + 0.1 


13.0 ±1.1 


11.8 ±2.2 


Mantz09 


Abell 2717 


0.0498 


2.6 + 0.1 




1.1 ±0.1 


Arnaud05 


Abell 3558 


0.048 


5.5 + 0.1 


7.7 ±0.6 


6.4 ± 1.0 


Mantz09 


Abell 1983 


0.0442 


2.2 + 0.1 




1.1 ±0.4 


Arnaud05 


MKW9 


0.0382 


2.4 + 0.2 




0.9 ± 0.2 


Arnaud05 
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Appendix B: L-z distribution of the cluster sample 
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Fig. B.l. Bolometric X-ray luminosity Lx of the clusters in- 
cluded in the combined cluster sample. 



Appendix C: Comparison of cluster properties for 
systems included in more than one subsample 
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Fig. C.l. Comparison of cluster properties for z < 0.3-systems 
included in more than one subsample. 

Several of the clusters in our combined sample are included 
in more than one source publication, either relying on com- 
pletely independent measurement data or using the same X-ray 
data reanalyzed according to the strategy that was chosen by 
the different authors. These measurements of cluster observables 
can be compared to each other after correcting for different anal- 
ysis schemes and provide a useful tool to test the applied homog- 
enization scheme. Furthermore, the comparison analysis gives 
an estimate of whether the error budget assumed in the differ- 
ent publications is realistic and reveals systematic differences 
between the results derived by the various studies. 

Fig. C.l shows the cluster properties of the low redshift 
(z < 0.3) overlap sample. In Fig. C.2, we present the deviations 



<-5a 



-4-3-2-10 1 2 3 4 



>5a 



^ 4 

i 3 



<-5a 



-4-3-2-101234 



>5a 



t5 7 

= 6 

O 5 

Z - I 



:-5o 




>5a 



ZhangOS 


Pratt09 


Arnaud05 


Mantz09 


Vikhlinin06 


Zhang07 


OHara07 



Fig. C.2. Distribution of deviations from the mean value for the 
z<0.3-clusters included in more than one subsample in units of 
the assumed error <x. Top panel: Mass deviations. Middle panel: 
Lx deviations. Bottom panel: ICM temperature deviations. 
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Fig. C.3. Comparison of cluster properties for z > 0.3-systems 
included in more than one subsample. 



of the individual measured values from the mean value in units 
of <x, and the error estimated in the individual publications. Note 
that the true values of the cluster observables are unknown. The 
comparison to the mean value therefore does not allow any state- 
ments about the reliability of the results derived in the differ- 
ent studies. Figure C.2 instead provides an insight into the sys- 
tematic differences between the results of different studies and 
whether the assumed error estimates are realistic. 

The derived spectroscopic temperatures agree well for most 
of the clusters. As visible in the bottom panel of Fig. C.2, the 
majority of the measured values deviate less than Icr from the 
mean value. In detail, 59% of the measurements lie within lcr 
and 82% within 2cr of the mean. Only 5% of the results deviate 
by more than 5cr. This indicates that the spectral fitting method 
generally leads to secure and consistent results, that the proba- 
bility of severe misestimations is low, and that the assumed er- 
ror budgets are likely to be realistic. Differences may result from 
different spectral extraction regions or different treatments of pa- 
rameters, such as the ICM metallicity or the background subtrac- 
tion process. However, these different measurement schemes do 
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Fig. C.4. Distribution of deviations from the mean value for the 
z> 0.3 -clusters included in more than one subsample in units of 
the assumed error cr. Top panel: Mass deviations. Middle panel: 
L x deviations. Bottom panel: ICM temperature deviations. 



not lead to completely incomparable data sets. The temperature 
differences between the subsamples are rather uncorrected and 
reveal no systematic trends between different studies. We note 
that for the samples of Zhang07, Zhang08, and Arnaud05 only 
core-excised temperatures were available. However, comparing 
those to the core-included temperatures given in other studies 
{e.g. Mantz09), the observed differences remain small. 

For the X-ray luminosity Lx, the situation is clearly different. 
As visible in the middle panel of Fig. C. 1 , most of the derived lu- 
minosities do not agree within the errors. Furthermore, the dif- 
ferences between the results of some studies clearly show sys- 
tematic trends. In terms of the deviations from the mean value, 
only 19% of the values lie within lcr and 29% within 2cr, while 
39% of the measurements show deviations of more than 5cr. The 
different samples exhibit systematic differences when compared 
to each other, especially for the Mantz09-Zhang08 overlap but 
to a lesser degree also for the common clusters of Zhang08 and 
OHara07. The reason for these deviations remains unclear since 
all known systematic differences, such as the definition of clus- 
ter radii and the different energy bands used, were corrected for. 
These deviations therefore imply that there are additional sys- 
tematic differences between the samples. However, for the cen- 
tral goal of this work, constraining the redshift evolution of scal- 
ing relations, this open question is of negligible importance be- 
cause systematic differences mostly occur for low-redshift sam- 
ples and the choice of sample from which multiply analyzed 
clusters are taken has no significant influence on the evolution 
results. In addition to systematic trends, even for samples that 
show no trends at all, the differences between the results consid- 
erably exceed the estimated errors. This indicates that the error 
estimations made by the different studies are too optimistic or 
that there are additional sources of measurement errors not in- 
cluded in the error budget. 

Similar but less significant systematic trends are also visible 
when comparing the results for cluster masses. As visible in the 
top panel of Fig. C.2, 51% of the results deviate by less than lcr 
from the mean value, while 89% lie within 2cr and no measure- 
ment shows deviations of more than 5<x. However, apart from 
these systematic trends the estimated errors for cluster mass 
seem realistic as most measurements deviate by less than lcr 
from the mean. The masses derived in Pratt09 based on the Yx- 



parameter and the Yx-M relation show no significant systematic 
difference from the hydrostatic mass estimates. However, ow- 
ing to the small overlap sample of five clusters, the comparison 
analysis provides no suitable tool to identify these differences. 

In Fig. C.3, the derived cluster properties of the z > 0.3 over- 
lap sample are compared. The typical observational errors for 
distant systems are larger, although within these errors the re- 
sults are in closer agreement than for the local overlap sample. 

The deviations from the mean temperature plotted in the bot- 
tom panel of Fig. C.4 lie below the estimated errors for most 
systems. In detail, 72% of the results deviate by less than lcr 
and 90% by less than 2cr from the mean value, while no devi- 
ations of more than 5<x occur. As for the low-redshift clusters, 
the measured temperatures show no systematic trends for single 
subsamples, i.e. the spectroscopic fitting procedure also seems 
reliable for distant clusters and there appears to be no major sys- 
tematic effects that have to be corrected. Furthermore, according 
to the mostly small deviations in units of the assumed error, the 
estimated error budget is likely to be realistic. 

Luminosities agree on average more strongly for the high- 
z clusters than for the local sample, for instance the Ettori04 
and OHara07 results are consistent for 12 of the 18 clusters in 
common (see middle panel of Fig. C.3). In contrast to the lo- 
cal overlap sample, no significant systematic trends between the 
different studies are visible. The deviations from the mean value 
plotted in the middle panel of Fig. C.4 are smaller than lcr for 
53% of the results and below 2cr for 62% of the measurements. 
We have found that 12% of the results deviate by more than 5cr. 
The distribution of deviations implies that the error budget might 
have been previously underestimated by the different studies, al- 
though by no means as significantly as for the local systems. 

The masses derived by Ettori04, Maughan06, and Kotov05 
plotted in the top panel of Fig. C.3 are consistent within the er- 
rors for all shared clusters, all results deviate by less then lcr 
from the mean value. The estimated errors are therefore likely 
realistic. However, the small size of the overlap sample of only 
four clusters does not allow us to peform a robust analysis of the 
systematic differences between the different studies. 

Appendix D: Local scaling relations for the 
combined cluster sample 

Fig. D.l, D.2, and D.3 show the z < 0.3-clusters included in the 
combined cluster samples and the local scaling relations fitted to 
this sample using different BCES fitting schemes in comparison 
to the relations derived by Pratt09 which were adopted for the 
evolution study in our work. 
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Fig. D.l. Local cluster sample: M-T relation. The red line shows 
the BCES(T|M) best-fit relation for the combined cluster sample, 
and the grey lines the BCES(M|T) and BCES orthogonal rela- 
tions. The blue line shows the Pratt09 relation (see Sect. 3.1). 
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Fig. D.3. Local cluster sample: M-Lx relation. The red line 
shows the BCES(L|M) best fit relation from the combined clus- 
ter sample, the grey lines the BCES(M|L) and BCES orthogonal 
relations. The blue line shows the Pratt09 relation (see Sect. 3.1). 
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Fig. D.2. Local cluster sample: L^-T relation. The red line 
shows the BCES(L|T) best-fit relation for the combined cluster 
sample, and the grey lines the BCES(T|L) and BCES orthogonal 
relations. The blue line shows the Pratt09 relation (see Sect. 3.1). 



